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ABSTRACT 


The  ARTY/GWC  method  for  providing  meteorological  messages  for  artil- 
lery firings  was  tested  in  this  study.  In  this  method.  Global  Weather 
Central  (GWC)  18-hour  prognostic  fields  are  updated  continuously  using 
Army  Corps  upper-air  soundings.  This  establishes  a prognostic  data  base 
from  which  meteorological  messages  may  be  extracted  when  required.  The 
testing  of  the  method  consisted  of  comparing  actual  and  simulated  cannon 
firings,  using  data  from  the  artillery  field  experiments  that  were  carried 
out  at  White  Sands  Missile  Range  in  November  and  December  of  1974.  The 
accuracies  of  the  simulated  artillery  firings  using  the  ARTY/GWC  method 
were  not  significantly  better  than  those  obtained  using  established  tech- 
niques. However,  improved  accuracies  with  this  method  should  be  possible 
when  prognostic  products  become  available  that  are  generated  by  numerical 
models  designed  to  treat  regions  of  limited  size  in  more  detail,  particu- 
larly with  regard  to  the  boundary  layer  and  mountainous  terrain. 
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I INTRODUCTION 


The  success  of  artillery  cannon  firings  depends  significantly  on 
the  accuracy  of  meteorological  (MET)  data.  For  the  purpose  of  improving 
on  ballistic  MET  messages  and  artillery  accuracies,  the  U.S.  Army  carried 
out  the  Prototype  Artillery  Subsystem  (PASS)  field  experiments  in  late 
1974  at  White  Sands  Missile  Range  (Barnett,  1976).  These  experiments 
consisted  of  Howitzer  firings  concurrent  with  upper-air  soundings  at 
ten  MET  stations  (Figure  1). 
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NAUTICAL  MILES 

FIGURE  1 UPPER-AIR  SOUNDING  STATIONS  AND  ARTILLERY  FIRING-IMPACT 
EMPLACEMENTS  — PASS  FIELD  EXPERIMENTS 
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The  PASS  data  were  used  recently  by  Blanco  and  Traylor  (1976)  to 
compare  different  methods  of  providing  artillery  MET  messages.  They 
found  that  methods  based  on  multistation  analyses  did  not  provide  signi- 
ficantly better  input  for  the  artillery  firings  than  did  the  standard 
method  of  using  upper-air  observations  of  a single  dedicated  station. 
However,  the  larger  firing  errors  were  during  days  when  the  wind  changed 
significantly  over  time  periods  of  one  to  three  hours.  This  indicates 
that  accurate  forecasts  of  the  changing  meteorological  conditions  should 
lead  to  significantly  improved  artillery  accuracy.  The  ARTY/GWC  method, 
tested  in  this  study,  is  an  attempt  to  provide  such  forecast  MET  mes- 
sages by  utilizing  the  U.S.  Air  Force  GWC  prognostic  products. 
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II  THE  ARTY/GWC  METHOD 

The  ARTY/GWC  method  was  developed  and  implemented  in  the  form  of 
computer  programs  in  a previous  study  (Mancuso  and  Hadfield,  1976). 

Flow  diagrams  of  the  various  computer  programs  used  in  the  present  ' 

study  are  shown  in  Figures  2 and  3.  The  GWC  prognostic  products  and  J 
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FIGURE  2 FLOW  DIAGRAMS  FOR  THE  DATA  PROCESSING  PROGRAMS  — 
GWCDAT  and  ARMDAT 
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U.S.  Army  upper-air  data  for  the  project  PASS  period  were  first  processed 
using  the  computer  programs  GWCDAT  and  ARMDAT  (Figure  2).  These  programs 
converted  vhe  data  into  a form  suitable  for  the  ARTY/GWC  programs.  The 
ARMDAT  routine  also  decodes,  identifies,  selects,  and  orders  (chronologi- 
cally) the  various  types  of  meteorological  sounding  data  that  were  col- 
lected during  project  PASS. 

The  computer  programs  of  the  ARTY/GWC  analysis  method  consist  of 
the  Prognostic  Data  Reanalysis  Routine  (PDRR)  and  the  Artillery  Applica- 
tions Routine  (AAR)  as  illustrated  in  Figure  3,  The  PDRR  uses  U.S.  Army 
weather  data  to  update  the  GWC  prognostic  fields  at  the  mandatory  levels. 
The  GWC  analysis  and  prognostic  data  are  used  initially  to  determine  the 
speed  and  direction  of  the  weather  pattern  movements  by  using  a pattern- 
tracking technique  developed  by  Wolf  (Mancuso  and  Wolf,  1974).  The  GWC 
data  are  then  interpolated  onto  a minigrid  that  provides  the  basic 
structure  for  performing  the  remaining  computations.  The  Army  upper-air 
sounding  data  are  used  to  update  the  GWC  prognostic  fields.  The  Army 
observations  are  grouped  within  40-minute  time  intervals,  and  are  treated 
as  sets  of  data  occurring  at  the  midpoint  of  the  time  intervals.  The  first 
set  of  observations  is  used  to  initially  update  the  GWC  prognostic  fields. 
This  updating  consists  of  obtaining  new  grid-point  values  based  on  both 
the  Army  upper-air  observations  and  GWC  grid-point  values.  The  new 
grid-point  values  are  computed  by  a least-squares  fitting  of  a first- 
degree  polynomial  to  nearby  data,  using  an  elliptical  distance  weighting 
(Mancuso  and  Wolf,  1974).  Difference  fields  between  the  original  and 
updated  GWC  prognostic  fields  are  then  computed.  These  difference  fields 
are  moved  with  the  speed  and  direction  of  the  GWC  patterns,  and  are  used 
to  update  the  GWC  prognostic  fields  at  later  times.  Second  and  subsequent 
sets  of  Army  observations  are  similarly  used  to  further  update  the  most 
recently  updated  GWC  fields. 

The  output  of  the  PDRR  program  (updated  GWC  fields)  is  used  as  input 
by  the  AAR  program.  In  the  AAR  computer  program,  an  updated  GWC  sounding 
profile  is  constructed  for  the  specified  time  and  location  of  the  cannon 
firing.  A guide  profile  is  also  used  to  aid  in  interpolating  between 
mandatory  levels.  The  guide  profile  is  obtained  either  from: 


4 


) 


5 


• A single  dedicated  station 

• The  best  available  station,  or  by 

• A multistation  analysis. 

The  latter  of  these  is  based  on  a first-degree  polynomial  fit  of  distance- 
weighted  data  as  in  the  updating  analysis.  However,  in  the  case  where 
only  the  simulated  Corps  stations  are  used  (Figure  1),  then  the  multi- 
station analysis  is  reduced  to  a simpler  distance-weighted  averaging 
of  the  data.  The  latter  is  similar  to  the  distance-  and  time-weighted 
algorithm  of  Blanco  and  Traylor  (1976).  The  AAR  program  is  also  used 
to  convert  the  final  updated  GWC  profile  into  a computer  MET  message 
format  (U.S.  Army,  1970).  The  output  of  the  AAR  program  (computer  MET 
message)  is  used  as  input  by  the  Ballistics  Research  Laboratory's  Bal- 
listic Trajectory  Simulation  (BRLBTS)  program. 


Various  improvements  and  corrections  to  the  previously  developed 
PDRR  and  AAR  programs  (Mancuso  and  Hadfield,  1976)  have  been  made  in 
this  study,  in  particular  to  give  an  optimum  analysis  and  editing  of 
the  data.  Listings  for  the  program  routines  that  have  been  modified 
are  given  in  Appendices  A and  B.  The  unlisted  routines  in  addition 
to  the  procedures  for  reading  in  program  control  parameters  are  identical 
to  those  given  previously. 
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Ill  PASS  AND  GWC  DATA 


The  PASS  upper-air  measurements  and  cannon  firings  were  made  at  the 
times  shown  in  Figure  4.  The  PASS  data  for  December  5th  and  7th  could 
not  be  used,  since  the  GWC  prognostic  data  were  available  only  through 
1800  GMT  of  December  5th.  Also,  the  initial  firings  on  November  12th 
and  14th  were  excluded  so  that  the  simulations  could  be  based  entirely 
on  the  GWC  prognostic  products  that  were  derived  from  1200  GMT  analyses. 
Thus,  a total  of  61  firing  cases  were  used  in  this  study. 

The  PASS  upper-air  measurements  were  made  at  the  ten  different 
stations  shown  in  Figure  1,  with  five  of  these  sites  being  used  to  simu- 
late an  operational  Corps  MET  network.  T9  radar  wind  measurements  were 
made  every  hour,  while  AN/GMD-IB  rawinsonde  measurements  were  made  every 
two  hours;  one-half  hour  separated  the  MET  observations  made  at  Corps 
sites  from  those  made  at  non-Corps  sites.  One  type  of  data  input  that 
was  used  in  this  study,  and  referred  to  as  GMD,  was  based  mainly  on  the 
rawinsonde  data.  However,  each  rawinsonde  observation  was  updated  at 
the  following  hour  by  replacing  the  rawinsonde  winds  with  the  more  recent 
T9  radar  winds.  A second  type  of  data  input  that  was  used,  and  referred 
to  as  GMD/T9,  consisted  of  winds  based  entirely  on  the  T9  radar  wind 
measurements.  However,  this  latter  type  of  input  was  inferior,  because 
the  T9  winds  were  less  complete  and  not  as  reliable  as  the  standard 
rawinsonde  winds. 

The  PASS  upper-air  measurements  had  been  converted  in  the  field  into 
both  computer  MET  messages  and  Air  Weather  Service  (AWS)  MET  messages. 
Computer  MET  messages  were  prepared  using  both  the  GMD  and  GMD/T9  type 
data.  However,  the  AWS  MET  messages  were  prepared  using  only  the  GMD/T9 
type  data. 
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IV  SIMULATION  APPROACH 


The  approach  used  in  this  study  to  evaluate  the  ARTY/GWC  method 
was  based  on  the  Ballistics  Trajectory  Simulation  program  developed  by 
the  U.S.  Army  Ballistics  Research  Laboratory.  With  this  program  the 
point  of  impact  is  computed  in  terms  of  the  range  (R)  and  a cross  com- 
ponent or  deflection  (D).  The  difference  between  a simulated  point  of 
impact  (Rg,  Dg)  and  an  actual  point  of  impact  (Rg,  D^)  is  caused  by 
errors  associated  with  the: 


• Interpolation  or  extrapolation  of  meteorological  data  in  both 
time  and  space 

• Meteorological  measurements 

• Location  of  the  actual  impact  (Rg,  Dg) 

• Quadrant  elevation  and  azimuth  angles 

• Muzzle  velocity 

• Projectile  characteristics  (weight,  shape,  etc.). 

Measured  values  for  these  last  four  elements  are  actually  based  on  mean 
or  representative  values  for  a series  of  cannon  shots  that  were  made 
within  a relatively  short  time  interval  about  a given  firing  time. 

In  order  to  evaluate  the  ARTY/GWC  method,  cannon  firing  simulations 
were  made  for  the  61  cases  using  MET  messages  from  the  ARTY/GWC  method 
and  from  more  basic  methods.  The  error  variability  of  any  given  set  of 
simulations  was  described  by  the  mean  and  standard  deviation  of  the 
quantities : 


AR  = Rg  - Rg  , 
AD  = Dg  - Ds  . 


(1) 


The  standard  deviations  of  these  quantities  and  and  a probable 

error  (PE)  are  used  primarily  for  comparing  methods,  assuming  that  the 
mean  deviations  either  are  due  to  other  sources  of  errors  or  would  ap- 
proach zero  if  the  sample  size  was  increased.  It  is  Important,  however. 
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to  consider  the  mean  departures  if  they  differ  significantly.  Error 
ellipses  were  constructed  for  comparison  purposes  using  standard  statis- 
tical theories  for  elliptical  frequency  distributions  (Brooks  and  Carru- 
thers,  1953).  The  probable  error  (PE)  as  defined  in  this  report  is  the 
radius  of  a circle  about  the  mean  which  contains  50%  of  the  impact  points. 
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V SIMULATION  RESULTS 


A.  GMD  Data 

The  initial  simulations  of  cannon  firings  were  made  using  a multi- 
station analysis  of  all  available  GMD  type  data  within  ±20  minutes  of 
the  firing  times.  This  is  equivalent  to  having  instant  data,  and  is  not 
realistic  for  Army  operations  with  present  measurement  systems.  However, 
such  simulations  demonstrate  the  accuracy  that  would  be  possible  with 
accurate  meteorological  forecasting.  The  multistation  analysis  was 
based  on  a first-degree  polynomial  fit  of  surrounding  data  weighted 
inversely  with  distance,  and  is  the  same  as  the  third  guide  profile  re- 
ferred to  in  Section  II.  The  results  for  these  initial  simulations  are 
shown  in  Figure  5a,  rtert  to  those  that  were  obtained  by  using  a standard 
atmosphere  with  zero  winds  (Figure  5b).  These  figures  show  the  differences 
between  the  actual  and  simulated  impact  locations  (AR,AD),  computed  ac- 
cording to  eqs  (1).  Since  the  location  of  the  simulated  impact  may  be 
thought  of  as  being  the  intended  or  desired  target,  the  figures  are 
similar  to  a typical  shooting  target.  As  discussed  previously,  errors 
from  several  sources  cause  the  differences  or  target  misses  that  are  dis- 
played in  Figures  5a  and  b.  The  simulation  results  that  were  based  on  the 
multistation  analysis  (Figure  5a)  have  a probable  error  (PE)  of  39m,  while 
the  results  for  the  standard  atmosphere  (Figure  5b)  have  a PE  value  of 
108m.  The  latter  also  has  a significant  mean  error  that  is  not  reflected 
in  the  PE  value.  These  two  results  illustrate  the  importance  of  meteorology 
in  artillery  ballistics. 

Four  different  tjrpes  of  simulation  results  are  shown  in  Figure  6, 
again  using  data  ±20  minutes  of  the  firing  times.  Figure  6a  repeats 
the  firing  simulation  results  for  the  multistation  analysis  that  were 
shown  in  Figure  6a  (using  an  expanded  scale),  while  Figure  6b  shows 
simulation  results  that  were  based  on  a simple  average  of  meteorological 
profiles  of  the  five  closest  stations.  This  latter  type  of  analysis 
gave  a firing  error  distribution  (Figure  6b)  very  similar  to  that  obtained 
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with  the  more  delineating  type  of  analysis  (Figure  6a).  This  suggests 
that  the  space  interpolation  errors  have  become  relatively  unimportant, 
and  that  the  scatter  shown  in  both  Figures  6a  and  b is  principally  due 
either  to  rawinsonde  measurement  errors  or  to  non-meteorological  errors. 

The  results  shown  in  Figures  6c  and  d are  for  cannon  firing  simula- 
tions that  were  based  on  MET  measurements  from  the  best  available*  sta- 
tion and  a dedicated  station  (station  1,  Fig.  1).  The  probable  errors 
(PE  “ 44  and  47m)  are  slightly  larger  than  those  of  Figures  6a  and  b.  The 
results  given  in  Figure  6 are  consistent  with  those  reported  by  Blanco 
and  Traylor  (1976)  that  were  based  on  only  rawinsonde  data. 

B,  GMD  and  GMD/T9  Data 

The  previous  results  were  all  based  on  the  GMD  type  of  data  that  were 
described  in  Section  III.  When  the  simulations  were  made  using  a multi- 
station analysis  of  the  GMD/T9  type  data,  the  cannon  firing  errors  in- 
creased by  about  20%,  as  illustrated  in  Figure  7.  Figures  7a  and  c show 
results  based  on  data  ±20  minutes  of  firings,  while  Figures  7b  and  d 
show  more  realistic  results  based  on  data  40  to  80  minutes  prior  to  the 
firing  time. 

C.  GWe  Data 

Computer  simulations  of  the  61  cannon  firings  were  also  made  using 
MET  messages  that  were  derived  completely  from  the  GWC  prognostic  fields. 
The  results  for  these  simulations  were  included  in  Figure  5 (c  and  d). 
Figure  5c  gives  the  results  that  were  obtained  based  on  the  GWC  products 
in  their  original  form.  The  probable  error  (106m)  is  comparable  to  that 
obtained  using  the  standard  atmosphere;  however,  the  mean  error  is  sig- 
nificantly lower.  Since  the  GWC  numerical  model  was  not  designed  for 
forecasting  mesoscale  features  (particularly  over  complex  terrain),  it 
was  believed  that  improved  results  might  be  obtained  by  replacing  the 


The  best  available  is  defined  as  the  observation  with  the  lowest  6, 
where  6 = (time  in  minutes  between  observation  and  firing)  plus 
(distance  in  kilometers  from  cannon  to  observation,  divided  by  2). 


aiaa 


FIGURE  7 DIFFERENCES  BETWEEN  ACTUAL  AND  SIMULATED  IMPACT  LOCATIONS, 
USING  A MULTISTATION  ANALYSIS  OF  GMD  TYPE  DATA  (a  and  b)  and 
GMD/T9  TYPE  DATA  (c  and  dl 
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wind  fields  at  the  surface  (<^80  mb)  and  at  850  mb  with  ones  that  were 
estimated  from  the  700  mb  wind  field.  To  account  in  an  approximate 
fashion  for  surface  friction  effect  on  wind  speed  (S)  and  direction  (6), 
these  estimations  were  made  as  follows: 


850  ” 2/3  Syoo* 

®850  " ®700  ’ 

sfc  " ®700  ’ 

®sfc  “ ®700  ' 

Figure  5d  shows  the  results  that  were  based  on  GWC  fields  modified  in 
this  manner.  Since  a slightly  lower  probable  error  (99m)  was  obtained 
using  these  modified  GWC  prognostic  products,  they  were  used  with  the 
ARTY/GWC  method. 

D.  GWC  and  GMD/T9  Data 

Computer  simulations  of  the  cannon  firings  were  made  based  on 
meteorological  input  as  provided  by  the  ARTY/GWC  method.  Since  the 
PDRR  program  of  the  ARTY/GWC  method  was  programmed  to  use  AWS  MET  mes- 
sages, the  simulations  were  made  using  the  GMD/T9  type  data: 

• At  the  time  (±20  min)  of  firing  (61  cannon  firings) 

• One  hour  (40  to  80  min)  prior  to  firing  (61  cannon  firings),  and 

• Two  hours  (100  to  140  min)  prior  to  firing  (48  cannon  firings). 

The  results  for  the  first  of  these  are  not  illustrated  since  they  were 
the  same  (except  for  minor  differences)  as  those  shown  in  Figure  7c; 
the  ARTY/GWC  and  multistation  methods  are  essentially  equivalent  when 
there  is  no  time  difference  (At)  between  the  MET  observations  and  cannon 
firings.  The  simulation  results  that  were  obtained  using  MET  data  both 
1 and  2 hours  prior  to  firings  (At  = -1  and  -2hr)  are  shown  in  Figures 
8a  and  b.  For  comparison  purposes,  simulation  results  are  also  shown 
for  a multistation  analysis  of  the  same  data  (Figures  8c  and  d).  The 
probable  errors  associated  with  the  ARTY/GWC  method  (PE  = 55  and  57m) 
were  sligjitly  larger  than  those  associated  with  the  multistation  method 
(PE  * 52  and  54m). 
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F IGURE  8 DIFFERENCES  BETWEEN  ACTUAL  AND  SIMULATED  IMPACT  LOCATIONS, 
USING  GMD/T9  TYPE  DATA 


In  the  above  experiments,  all  of  the  MET  stations  in  Figure  1 were 
utilized.  When  the  more  realistic  set  of  simulated  Corps  stations  de- 
signated in  Figure  1 was  used,  then  the  results  shown  in  Figure  9 were 
obtained.  In  these  simulation  experiments,  the  ARTY/GWC  gave  comparable 
results  to  those  of  the  multistation  method  (54  and  56m  versus  53  and 
57m),  but  the  mean  errors  for  the  ARTY/GWC  method  were  slightly  larger. 

In  the  results  shown  in  both  Figures  9 and  8,  the  ARTY/GWC  method  actually 
provided  superior  temperature  MET  messages  which  resulted  in  lower  range 
errors  but  larger  deflection  errors,  since  temperature  affects  only  the 
range. 

E.  GWC  and  GMD  Data 

The  GWC  prognostic  products  influence  the  final  artillery  MET  mes- 
sage in  the  following  three  ways: 

• The  GWC  values  are  included  in  the  update  analyses,  but  with  a 
weighting  that  becomes  significant  only  at  large  distances 
from  the  Corps  data. 

• The  updated  GWC  fields  are  displaced  in  time  using  the  motions 
of  the  GWC  patterns. 

• The  time  changes  in  the  GWC  fields  are  added  directly  onto 
the  displaced  updated  fields  at  later  times. 

Inspection  of  the  individual  cases  showed  that  the  ARTY/GWC  approach 
failed  when  there  were  large  time  changes  in  the  wind  fields  predicted 
by  GWC.  This  occurred  noticeably  on  December  2nd,  when  the  GWC  prognostic 
fields  showed  suspiciously  large  changes  over  a three  hour  period,  that 
were  inconsistent  with  the  Corps  rawinsonde  data. 

Therefore,  final  computer  simulations  were  made  in  which  the  third 
of  GWC  Influences  mentioned  above  was  excluded  from  the  ARTY/GWC  method. 

In  these  final  simulations,  the  GMD  data  were  used  rather  than  the  GMD/T9, 
since  the  errors  associated  with  the  latter  tend  to  obscure  the  results. 
This  required  a special  program  for  converting  computer  MET  messages 
into  AWS  MET  messages.  The  firing  errors  that  were  associated  with  using 
data  both  1 and  2 hours  prior  to  the  firings  (At  = -Ihr  and  -2hr)  are 
shown  in  Figures  10a  to  d.  The  ARTY/GWC  forecast  method  gave  only  slightly 
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FIGURE  9 DIFFERENCES  BETWEEN  ACTUAL  AND  SIMULATED  IMPACT  LOCATIONS,  USING 
GMD/T9  DATA,  RESTRICTED  TO  CORPS  STATIONS  (See  Figure  I) 
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smaller  probable  errors  than  did  the  multistation  persistence  method 
(47  and  52m  versus  48  and  53m). 


The  AR  and  AD  values  for  individual  cases  are  listed  in  Appendix  C 
for  the  results  shown  in  Figures  5,  6,  and  10. 
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VI  CONCLUSIONS  AND  RECOMMENDATIONS 

This  study  showed  that  a significant  reduction  (25%  to  50%)  in  artil- 
lery firing  errors  is  possible,  if  wind  and  temperature  profiles  for  the 
time  and  location  of  the  artillery  firing  can  be  accurately  determined, 
either  by  direct  measurements  or  by  accurate  forecasting.  The  ARTY/GWC 
method  for  providing  forecast  artillery  MET  messages,  that  was  tested, 
provided  results  statistically  equivalent  to  simpler  persistence  methods. 
This  indicates  that  the  GWC  prognostic  products  are  not  yet  sufficiently 
accurate  and  detailed  for  use  in  artillery  ballistics.  If  more  advanced 
GWC  models  become  available,  which  are  designed  to  treat  limited  and  re- 
locatable window  regions  on  a mesoscale  basis,  then  the  ARTY/GWC  or  similar 
method  could  prove  to  be  highly  effective. 

In  this  study,  no  attempt  was  made  to  include  the  influence  of  com- 
plex terrain  in  the  ARTY/GWC  updating  analyses.  In  actual  military  oper- 
ations, this  would  frequently  be  of  greater  importance  than  it  is  in  the 
White  Sands  areas.  Theoretical  and  diagnostic  techniques  such  as  those 
that  have  been  developed  by  Estoque  and  Bhumralkar  (see  Estoque  et  al., 
1976),  can  be  used  to  investigate  the  influence  that  terrain  has  upon 
the  meteorological  fields,  with  regard  both  to  roughness  and  surface 
heating.  Research  along  this  line  is  currently  being  carried  out  at 
SRI  International. 
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Appendix  A 


PROGRAM 

2 


LISTING  OF  MODIFIED  PDRR  PROGRAM  < 

I 

PORR< INPUT. OUTPUT, TAPEl .TAPE2, TAPE3. TAPE4 .TAPESs INPUT . 

TAPE9, PUNCH) 


C THIS  PROGRAM  IS  THE  PROGNOSTIC  DATA  REANALYSIS  ROUTINE  (PDRR).  IT  IS 
C USED  FOR  REANALYZING  OR  UPDATING  THE  GLOBAL  HEATHER  CENTER  I GHC ) 

C PROGNOSTIC  DATA  USING  THE  MOST  RECENTLY  AVAILABLE  ARMY  UPPER  AIR 
C OBSERVATIONS. 

C AUG  1977  version:  - CONTAINS  IMPROVEMENTS  SUCH  AS  TEMPERATURE  EDITING. 
C OPTIMIZED  PARAMETER  VALUES.  AND  MORE  ACCURATE  COMPUTATION  OF  CMC 
C PATTERN  MOTIONS. 

C 

C DIMENSION  AND  COMMON  STATMENTS 
C 

INTEGER  DATA! 12000) 

DIMENSION  U( 1 00  ) .V( 100) .HI  100) .TIIOO ) . VORI 100) .01  VC  100) .BALC 100 ) 
DIMENSION  use  2S).VS(2S).HS(25) , rSC  25) . XS(25) . YSC  25) .ESC  25) 
DIMENSION  XC 10)  ,YC 10) . OS X ( 1 00 ) . OSY I 1 00 ) .LD ATE ( 5 . 25 ) . I MDC 12) 
DIMENSION  CMC  10) .TMC 10) . TRUC 1 00  ) .TRVC 1 00) .XRUC 1 00)  .XRVC 100 ) 
DIMENSION  OAT1C400 ) .GATl C 100 ) .GAT21 1 00 ) . V AT C 25 ) . MA TC 2 5 ) .GASC 100.4) 
DIMENSION  URC 100) . VRC 100) .HRC 100 ).TRC 100) 

COMMON/CIS/  ISC  1000) 

COMMON/CCHK/  CR T . S 1 M .G I M . 1 WNO 
COMMON/CGD/  M9,N9. I 9 . MB , N8, YB . XB . DO 
COMMON/CPTS/  KS.Wl .C2.RMAX.KSSS. IDS.KSW.ALPH 

EQUIVALENCE  CU.GASC 1 . 1 ) ) . C V.GASC 1 .2) ). C H.GASC 1.3)) . IT. GASC 1.4)). 

2 1UR.OAT1C1)).CVR.OAT1I101)).1HR.OAT1C201> ). 

3 C TR.  OATH  301  ) ) . C VOR.OSX  ) . C DI V . DSY  ) 

C 

C SET  BASIC  CONSTANTS  I INITIALLY  SET  CORE  TO  ZERO) 

C 

DATA  XNIL.TMAX. TMIN/-999.9.310.0.230.0/ 

DATA  I8CK,  IDIF.IGMT  /O. 20. -420/ 

DATA  10.11.12/0.1.2/ 

DATA  NVARP.NSIZE  /3.5/ 

DATA  CRT.SIM.GIMS.IWNO/  0.3.0.10.0.05.1/ 

DATA  N1  .C2 .RMAX .K S . KSS5 , KSH . I OS . ALPH/0 . 0 1 . 0 . 0 05. 10.0.10.5.4.1.2.0/ 
DATA  I20.L1 .L4.L1 00.K4.N5.N2S  / 1 6 . 1 . 4 . 1 00 . 4 . 5 . 25/ 

DATA  XP.YP.ROT. OTS. ACR/24.0.26. 0.80. 0. 1 .0.0.0174533/ 

DATA  I MO/0. 31 .59.90.120. 151 . 181.212.243.273.304.334/ 

1 FORMAT! lOX. 71  10  ) 

2 FORMAT! 10X,2I 10,5F10.2) 

3 FORMAT! lOX. I 1 0) 

4 FORMAT! 14X . 3I2.6X.2  12) 

5 FORMAT! lOX. I 10) 

6 FORMATC10X.3F10.2) 

7 FORMATC2X.6I 1 0.2E1S.3) 

8 F0RMATC2X.4E1S.3) 

9 FORMATC2X.6H0ELETE.I10.4F7.2) 

10  FORMAT! 10X.5F10. 2) 

C 

C READ  IN  CONTROL  PARAMETERS 
C 

READ  !5.1)  ICHK. I ANA. IDI V. IBAL. IFOR. IPRINT. I TAPE 
READ  !S.2)N9.M9.XB.YB.D0.XC.VC 
READ  IS. 5)  JT 

READ  C5.6)  eXS!  J).rS!  Jl.ESCJ)  .3.^1.  JT) 
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I 
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) 
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I 

I 

I 
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I 

f 


PRINT  OUT  CONTROL  PARAMETERS 


PRINT  9000 
PRINT  9100 
PRINT  9000 
PRINT  9050 

PRINT  9001 . ICHK.I ANA. 101 V.IUAL. IFOR* I PRINT* 1 TAPE 
PRINT  9002«N9*M9*XB*YB*00*XC*VC 
PRINT  9005.JT 
PRINT  9056 

PRINT  9006* (J  *XS(J ) *VS(J 1 *ES4Jl * J=1 * JT ) 

PERFORM  BASIC  COMPUTATIONS 

IF  (JT.CT.25)  GO  TO  999 
IF  (M9.GT.10.QR.N9.GT.10)  GO  TO  999 
ACRI^l .0/ACR 
GlMsGIMS 
XI 00=OTS*60.0 
I9aM9*N9 
M8*M9-1 
N8=N9-1 
N49=I9 
N98=NA94-N49 
N47=N984N49 
K3=K4-1 
I0=L4«K3*N49 
DO  20  N=1*N9 
20  X<N)3XB400*IN-1 ) 

XAVE= ACR*(-X( N5)-R0T) 

SAVEsSINIXAVE ) 

CAVEaCOSI XAVE) 

DO  22  M=1*M9 
VI M)=Ya-00*I M-1 1 
ANG=ACR4Y(M I 
CM(M)=COSIANG) 

22  TM(M)zTAN< ANG) 

XNS=-X(NS) 

CGO-1 «04V(NS) /90.0 

CALL  MESHI I 0* JT * VS* XS* VS. US *HS* TS. V . X. V *U.H.T ) 

IH20=I 120-1 >460 
OO  23  1=1*100 
TRUI I 1=0.0 

23  TRV(1)=0.0 
1=0 

OO  24  M=1*M9 

TL A=TAN( (90.0-Y(M))*ACR*0.S>431.2042 

OO  24  N=1 *N9 

I=141 

XL0=ACR4 I -X  I N )-R0T ) 

OSXI I l=IXP-TLA«SIN<XLa)-XC)*2.041 .0 

24  OSYf I >=4VP4TLA4C0S(XL0>-YC)42.041.0 
RE AO  45*3)  IT 
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READ  (S.4)  ICLOATE<L*n<L=lf5>.l  = l f IT) 
PRINT  9003. IT 
PRINT  9054 

PRINT  9004.  KLOATEIL.  I ).C-I.5).  1 = 1 .IT) 

C 

C 

C LOOP  THROUGH  ALL  GWC  FORECASTING  TIMES 
C 

DO  50  IG=1.I20 
READ  IN  G«C  DATA 

CALL  CWCINIIG.IGWC. JOATE. JTIME.OATl. IMO) 
31  DO  45  L=L1.L4 
Nl=IL-l )*L100 
IP=< I IG-I )*L4*L-I )4I9*K3 
JP=0 

DO  42  K=1.K4 
IF  IK.EQ.K4)  1P=IP-N49 
DO  36  N=I.N25 
NI=N141 
OATN=OATItNI) 

«AT(N)=OATN 
IF  (K.GT.2)  GO  TO  35 
N3=N1 

IF  IK.E0.2)  N3=N1-N25 
OArN=OATt(N3) 

VAtN*OATl <N34N25) 

MAT(N)sOATN 

IF  (K.E0.2)  MATIN)=VArN 
IF  fK.EQ.2)  OATN-=VATN 

35  JP=JP4I 
GATI I JP)=CAS( JP.L) 

GASI JP,L)=OATN 
GAT2( JP)sOATN 

36  VAT|N)=DATN 
DO  40  1=1.19 
OLX=OSX( I ) 

OLV=DSV< I ) 

NX=INT(OLX) 

Mr=INT<OLY) 

NM=( MY-1 )•N5♦NX 
OLX=OLX-NX 
OLY=OLY-MY 

C 

C INTERPOLATE  TO  MINI  GRID 
C 

CALL  1NTPTINM.N5.VAT.0LX.0LY.BAT) 

IP=IP41 

1BAT=1NT(BAT*I0.0) 

IF  (K.E0.3)  IBAT=iaAT*I0000 
1F(K.NE.4)  GO  TO  40 
IFIDATAI IP) .LT.O)  tBAT=>lBAT 
1BAT=0ATA1 IP)4iaAT 
40  DATA! |P>=iaAT 


A-4 


1 

f 


no  non 


42  CONTINUE 


COMPUTE  TRENDS 

IF  ( IFOR.LT.l .0R.IG.LT.2)  60  TO  45 
LIGaf  IG-1  »♦L4♦^. 

NVARSsNVARP 

CALL  TREND! NVARSi NS IZE.TRUILIG) tTRyCLIG) .GATl .6AT2. IG ) 

TRUfLlG>-TRU(LIG)«CGO 

TRVIL1G)»-TRV|LIG)*CGD 

CALL  UVCONVITRUILIG) tTRVILIGI .XN5) 

TRU!LIG)=TRU(LIG)/00 

TRV!LIG)=TRV!LI6)/DO 

TRUV=TRU(LIGI *TRU(L IG)4TRV(LIG)«TRVfLIG> 

IF  (TRUV.LT. 16.0)  GO  TO  44 
TRFC-4.0ySQRTITRUV> 

TRU(LlG)sTRU{LIG)*TRFC 

TRV<LIG)aTRV(LIG)*TRFC 

44  CONTINUE 

45  CONTINUE 
50  CONTINUE 

ILT*I20*L4 
DO  55  LSL1.L4 
DO  55  IG-S2.I20 
LI6»( IG-I )4L44L 
SUM3;2.0 

USUM=TRU( L I G I 42 . 0 

VSUMsTRV(LIG)42.0 

DO  54  IGX=1.3 

XIGX^l .O/IOX 

LI0sl1G-I’IGX*L4 

IF  (LID.GT.ILT)  GO  TO  53 

SUM=SUMTX IGX 

USUM=USUM'»’TRU(LI0)*XI6X 

VSUM^VSUM^TRVILIOIVXIGX 

53  LI03LIG-IGX4L4 

IF  (LIO.LT.I)  GO  TO  54 
SUM^SUMTX IGX 
USUMsUSUM*TRUtL ID) *X IGX 
VSUMaVSUM4TRV(L ID) •XIGX 

54  CONTINUE 
SUM«t .0/SUM 
XRU(LI6>«USUM*SUM 
XRV(LIG)«VSUM*SUM 

55  CONTINUE 
PRINT  9025 

DO  56  L-L1.L4 
DO  56  I6S2.I20 
LIG=( IG-I )*L4*L 
TRUCLIGIsXRUtLlG) 

TRV(L1G>3XRV(LI6) 

PRINT  2. IG.L.TRU(LIG) .TRVILIG. 

56  CONTINUE 

END  OF  G«C  LOOP 
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PRINT  2.JP.IP 


LOOP  THROUGH  ARMY  RAWINSONOE  OBSERVATION  TIMES 

DO  100  IRsl.IT 

IF  llANA.LT.ll  GO  TO  101 

1*IR 

L02=LDATE12.I  I 

lOATEsLOATEtS.I  l-fOO*!  LDATEI  4 1 I >424  4<  LO  ATEI  3 • I )«'IM0(L02> 
2 «’365*LDATEfl  •11-366  M-IBCK 

PRINT  1»1GMC*I0ATE 
IF  I 10ATE.LT. IGWCl  GO  TO  999 
IF  (IOATE.GT.16MC41H20)  GO  TO  999 
1COMP3IGWC4X100 
OO  60  16*2,120 
IF  ( lOATE-lCOMPl  56.60.60 
56  FGl-< ICOMP-IOATEI40.1/XIOO 
FG2=0.1-FG1 
GO  TO  61 

60  1COMP3ICOMP4X100 
GO  TO  999 

61  1C=IG-1 
IGl3lG4i 
IL=I lG-1 >410 

PRINT  1 . ICOMP.IDATE. IG 
PRINT  6.FG1.FG2 

READ  IN  RAMINSONOE  DATA 

I=IR4IGO-l 

CALL  RAMINljr.I.  ES .UR • VR .HR . TR. IMO. lOATE. lOIF • IGMT 1 

LOOP  THROUGH  MANDATORY  LEVELS 

OO  95  L-L1.L4 
LJ=IL-1 > 4JT 

00  68  J-l.JT 
US(J)«UR<  J4LJ  1 
VSI J)=VR< J4LJ) 

HS(J»»HR( J4LJ 1 

68  TS(  J MTRI  J4LJ  1 

N1=IL4<L-1 I4N494K341 
N2=N14N49-1 

1 = 1 

DO  70  N=N1.N2 

U( I |30ATA|N|4FG140ArA(N4I0|4FG2 
VI 1 tsOATAIN^NAO  ) 4FG 1 40 AT A< N4N494 lO 1 4FG2 
IHl=tNT(DATA(N4N98  140.0001 1 
1H2=INT(OATA(N4N984I0140.0001 1 
1T1S0ATACN4N98  1-IH1410000 
I T2*0ATAIN4N984I0)- IH241 0000 
tTl=lABS( ITl ) 

IT2=IABS1 IT21 
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H(I)3  IH1*FG1 «IH2*FG2 
T<II»  IT1*FG14’IT2*FG2 

70  I=I*l 
PRINT  9000 
PRINT  9059.1. 

SUMTi^O.O 

NUMT=0 

DO  166  J^l.JT 

IF  (TSIJI.LT.TMIN  .OR .TS <J > .GT . TMAX  ) GO  TO  166 

SUMTsSUMT^TSf  Jl 

NUMT=NUMT-»1 

166  CONTINUE 

IF  (NUMT.LT.l)  GO  TO  169 

SUMT=SUMT/NUMT 

OO  168  Jsl.JT 

ABTS^ABSITSIJ  )-SUMTl 

IF  I ABTS.LT.8.0 1 GO  TO  168 

IF  ITSIJI.LT.O.OI  GO  TO  167 

PRINT  9.J.TS(J).HS(JI.USIJI.VS(J) 

167  TS<J»=XNIL 
HSI J)=XNIL 
US< JlsXNiL 
VS( JlsXNIL 

168  CONTINUE 

169  CONTINUE 

IF  ( ICHK.LE.O.OR.L.EQ.I 1 GO  TO  170 

CALL  CHECK! JT.YS.XS.VS.US.V.X.V.U) 

170  CONTINUE 
PRINT  9057 

PRINT  900  7.! J.HSf  J).TSIJ>.US!J).VS1J).J=1.JT> 

IF  ! IPRINT.lt. 1 ) GO  TO  71 
PRINT  9060 
PRINT  9062 

PRINT  9010.! H! I > . 1=1 . 19) 

PRINT  9064 

PRINT  9010. ! T! I ) . 1 = 1 .19) 

PRINT  9066 

PRINT  9010.!U!I>.1=1.I9) 

PRINT  9068 

PRINT  9010. ! V! I ) . 1 = 1 . 19) 

71  CONTINUE 

IF  !IANA)  95.95.73 
73  OO  173  1=1.19 
VOR! I )=0.0 
OIV! I )=0.0 
173  BAL!I>=0.0 

IF  IIOIV.GT.O) 

1 CALL  KID! 10. I 1 . 10 . VOR.OI V .BAL .U.V.CM.TM) 

CALL  MESHII2.JT.VS.XS.VS.US.HS.TS.Y.X.V.U.H.TI 
IF  IIOIV.LT.O)  GO  TO  75 

ADJUSTMENT  OF  WIND  FIELD  TO  THE  GMC  DIVERGENCE  FIELD 

CALL  K ID! I 1 . 10. 10* VOR. BAL .01 V.U. V.CM . TM) 

CALL  ALTERSI25.0.0.0.5.TM.  CM.  U.V.VOR.OIV) 


non 


75  IF  (tBAL*LT.l)  GO  TO  85 

COMPUTATION  OF  BALANCED  HEIGHT  FIELD 

CALL  KlOIIOtlO.II •VOR«OIV.BALtUtV.CM.TM) 

SUM=0*0 
DO  80  I=I . 19 
80  SUM=SUM'*’H(I) 

AVE=SUM/I9 

CALL  BALHGTI M9«N9>0D> !• 10*1 .2*  « 01 • A VE (CM . TM .H t BAL ) 

85  CONTINUE 

IF  < IPRINT.lt. 1 > GO  TO  86 
PRINT  9000 
PRINT  9070 
PRINT  9062 

PRINT  9010. (HI  I I . I=I . 19) 

PRINT  9064 

PRINT  90IO.(T(  I ) . I-I .19) 

PRINT  9066 

PRINT  9010. (U( I ). I=I . 19) 

PRINT  9068 

PRINT  9010. I V( I ) . I-l. 19) 

86  CONTINUE 

IF  (IFOR.LT.l)  GO  TO  95 
C 

C UPDATE  ALL  GMC  FORECAST  FIELDS  BT  ADVECTING  MITH  COMPUTED  TRENDS 
C 

DO  90  IH3IG.I20 

NI=( IH-1 )*ID4IL-l )4K3*N49+l 

N2=N14N49-1 

OT=OTS 

LlHs( IH-l )*L44L 
IF  (IH.GT.IGI)  GO  TO  88 
LIK=LIH4L4 
OT--FG2410.04DTS 
IF  (IH.EQ.IG)  GO  TO  187 
OT-  FG1410.04DTS 
DO  186  Isl.19 
U( I )=XRUf I ) 

VI  I )3XRVI  I ) 

HI  I ) = VOR( 1 ) 

186  TII)=:0IVI1) 

GO  TO  88 

187  1=1 

DO  87  N=N1.N2 

UI  1 )=0ATAIN)*FG1«’DATAIN4I0)  4FG2-UI  I ) 

VII  )=0ATAIN4N49)«FG140AT  AIN4N494'  1D)*FG2-V(  I ) 

IHI=INTIOATAI N4N9a  )*0.0001) 

IH2=INTIOATA(N4N98410I*0 .0001 ) 

ITI=OATA(N4N98  >-lHl*10000 
ITEsOATAIN'fNOS^  ID )•  1H2*  1 0000 

ITlsIABSI ITl ) 
lT2=IAaSI IT2) 

HII)=  IH1*FG1-»IH2*FG2-HI  I ) 

TII)=  IT14FG14IT24FG2-TI I ) 
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c 

c 

c 

c 

c 

c 


XRUi 1 )sU(  I ) 

XRV( I ) = V< 1 ) 

VORC I )=H( 1 ) 

Dive  I )=T( I ) 

87  I = I«-1 

88  CONTINUE 

CALL  AOVECIOT.TRUILIH) .TRV(LIH) .UtBAL.CMI 
CALL  AOVEC(OT*TRU(LIH) .TRV(LIH) .V.BALtCM) 
CALL  ADVECIDT  .TRU(LIH) .TRV(LIH) .H.BALtCMI 
CALL  AOVECIOT  tTRUILIH) .TRyiLIH) .T«BAL.CH) 
I=l 

DO  89  NSN1.N2 

DATA!  NtsOATAI  N)>U(n*10 
0ATAlN«'N49)=0ATAf  N-»N49)-V(  I )*10 
IHl=INT(OATAC N4N98  )*0.0001) 

IT1=DATA(N4N98  )-IHl*lOOOO 
ITl=IABS( ITl ) 

IH2-IH1-H( I )*10.0 
IT2=IT1-T(I )*10.0 
IF(  IH2.lt. 0)  IT2=-IT2 
OATA(N'*'N98I-IH2*IOOOO'«’IT2 

89  1=141 

90  CONTINUE 
95  CONTINUE 

END  OF  LEVEL  LOOP 


WRITE  OUT  UPDATAO  FIELDS  ON  TAPES 
101  CONTINUE 

IF  I 1TAPE.lt.  1>  GO  TO  99 
PRINT  I.  JOATE. JTIME, 120 
WRITE  <3)  JOATE. JT IME. 120 
DO  98  IH=I.I20 
DO  98  L=L1.L4 

N1=I IM-1 )4ID4(L-1 >*K3*N4941 
N2=N14N49-1 

98  WRITE  (3)  lOATA(N) .0ATA(N4N49).0ATA(N4N98) .N=N1.N2) 

99  CONTINUE 
100  CONTINUE 


END  OF  ARMY  LOOP 


END  FILE  3 


197 

READ  (11  IDUM 

IF  (EOF(ll)  200. 

197 

200 

CONTINUE 

GO  TO  1000 

999 

PRINT  9999 

1000 

CONTINUE 

9000 

FORMAT! IHl 1 

9001 

FORMAT(2X.8HCARD 

A 

.7I10//> 

9002 

F0RMAT(2X.aHCAR0 

B 

•2110.SF10.2//) 

9003 

FORMAT<2X.8HCARO 

C 

• 1 10//) 

\ 

{ 


9004  FORMATCSX.SI 1 0/ I 

9005  F0RMAT(/2X.8HCA«0  E ,110.//) 

9006  FORMATi lOX. IS.SX.3F10.2/ ) 

9007  FORMAT! 16Xt 12 .2X,4F10.2) 

9010  FORMAT<IP.10X.7F6*0> 

9025  FORMAT! IH  ) 

9050  FORMAT!30X,27HINPUT  DATA  FOR  PDRR  HOUTINE///1 

9054  FORMAT!2X.57HCAROS  0 YEAR  MONTH  DAY  HOUR 

2 MIN/> 

9056  FQRMAT(2X»aHCARDS  F . 9X • 20HST AT  I ON  t NFORMAT I ON// 1 2X , 1 OHNUMBEri 

2 ,30HLONGI TOTE  LATITUTE  ELEVATION/) 

9057  FORHAT!25X.t3HSrAriaN  DATA//  9X,33H  NUMBER  D VALUE  TEMPERAT 

2URE.  19H  U COMP  V COMP  ./) 

9059  FORMAT! 10X,22HRESULTS  FOR  LEVEL  L = •ll//> 

9060  FORMAT!20X,2 IHGWC  DATA  ON  MINI  GRID//) 

9062  FORMAT1/24X.8HO  VALUES/) 

9064  FORMAT(/23X.l 1 HTEMPERATURE/ ) 

9066  FORMAT1/26X.6HU  COMP/) 

9068  F0RMAT!/26X*6HV  COMP/) 

9070  FORMAT! 15X,29HUPOATEO  GWC  DATA  ON  MINI  GRID//) 

9100  F0RMATI///45X .34HPROGNOSTIC  DATA  REANALYSIS  ROUTINE.////// 

2///.  57X,5HUNI TS.///.53X,l IHSPEEO  - MPS . // , 53X , 1 5H0 IRECT I ON  - DE 
3G«//.  53X, ISHHEIGHT  - METERS. //, 53X , 1 9HTEMPERATURE  - DEG  K) 

9999  FORMAT! IOX.23HINCONSISTENT  INPUT  DATA) 

STOP 

END 


non  non 


SUBROUTINE  RAWIN( JT. I t E S • UR t VR . HR • TR . I MD • IDA TE. lO IF . IGMT > 

C 

C THIS  SUBROUTINE  READS  IN  THE  U S ARMY  RAMINSONOE  DATA. 

C AUG  1977  version:  - MODIFIED  TO  SELECT  DATA  TYPES  IT=4  OR  5 ONLY,  AND 
C TO  INTERPOLATE  AND  EXTRAPOLATE  WHEN  WIND  PROFILES  ARE  INCOMPLETE. 

C 

DIMENSION  UR(1  ) . VR( I ) .HR! 1 I ,TR( 1 ) ,ES( 1 > . IMOI 1 ) 

DIMENSION  OAT(S12).IS(2SI .HOZ( 17) ,PRS( 6) .STD (6) ,C0R(6) 

DIMENSION  IR(SO) 

DATA  ACR«CKM.CC.XNIL/0*0I74533.0*5148.29.29.-999.9/ 

DATA  JIO. J17. J24, J28. J51 .K4/10, 12, 24. 28. SI, 4/ 

DATA  STD/  0.0. 1457.0,3012. 0,5574.0,0. 0.0.0/ 

DATA  PRS/  870.0.850.0.700.0.500.0.300.0.200.0/ 

DATA  HOZ/0.0 .100.0 .350.0 .750.0. 1250.0. 1750.0,2250.0.2750.0,3250.0 
2 .3750.0.4250.0.4750.0.5500.0.6500.0.7500.0.8500.0.9500.0/ 

DATA  COR/-3.e.-3. 1 .-0.1. 0.5.0. 0.0.0/ 

SPDMX= 100.0 
ESST=1300.0 
DO  18  L= 1 . J T 

18  IS<L)=11 
JK=  JT4X4 
DO  19  L=1.JK 
1R1L)=0 
HRlDsXNlL 
TR(L)=XNIL 
URID^XNIL 

19  VRID^XNIL 
IFIL=0 

IF  (I .CT.l ) GG  TO  22 
READS  IN  DATA  SET 

20  READ  (9)  OAT 
IF  (E0F<9))124,21 

124  IFIL=IFIL«-l 

IF  (IFIL.GT.l)  GO  TO  32 
GO  TO  20 

LOOP  THROUGH  DATA  SET  IN  RECORD 

21  CONTINUE 
DO  30  J=1.J10 

22  JC=lJ-l)*J5141 
C 

C CHECKS  DATA  FOR  0 ATE , T I ME , AND  TYPE 
C 

IF  (DAT< JC)  .LT.0.0  ) GO  TO  30 
IP1  = INT<DAT( JC) *1 .OE-4  ) 

DATJ=OAT( JC)-IP1*1 .0E4 
IP2=INT<OATJ*l  .OE-2  ) 

IP3  =OAT J-IP2*1  .0E2 
I P4=  I NT  (OAT  (JO  I >40.01  ) 

1P5=DAT( JC41 )-IP4*100 

JOATE»IP5460*I  IP4>24  4(  IP34IM0(  IP2  ) ♦•365*  IPl- 366  ) ) - 1 GMT 
5 FORMAT(2X.SI 12) 
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IF  t JDATE.lt, IDATE-IDIF)  GO  TO  30 
IF  ( JDATE.GT  . IDATE+ IDIF)  GO  TO  32 
JS=INT(OAT( JC  + 2IF0.1  ) 

IT  =DAT( JCF2 >-JS*lO 
IF  ( IT.LT.4.0F.IT.GT.6>  GO  TO  30 
NOATE=INTtDAT( JC) ) 
NTIME=INT10AT€JCH I) 

ESJS^ESI JS) 


SELECTS  OUT  HEIGHTS  AND  TEMPERATURES  FCR  THE  MANDATORY  LEVELS 
K = I 

DO  26  JQ=A,J24,4 
JK=<K-1»4JT  ♦JS 
JCO=JC+JQ 

IF  ( K.EQ. 1 ) GO  TO  25 

IF  (K.GT.K4)  GO  TO  27 

IF  (OAT< JCO) .NE.PRSiK) > GO  TO  26 

25  IF  (OAT( JCQ-1  I.LT.O.O)  GO  TO  125 

IF  ( IT.GT  . ISC JS ) .AND.TRI JK ) .GT.0.0)  GO  TO  125 
HR(  JK  )-=OATC  JCQ-  I )-STD(K> 

TRC JK»=OAT( JCQ+l ) ♦273. 16 
IF  (K.NE.l)  GO  TO  125 

PRSl=(<ESJS-ESST)/(TRC JK  > 4CC > ♦ 1 . 0 ) *0 AT ( JCO ) 

HfiC JK)=-(DAT( JCQ)-  PRSl  ) ATR ( JK I *CC/0 AT C J CQ ) 

125  K=K*l 

26  CONTINUE 

27  KT=K-1 


INTERPOLATES  ZONAL  4 1 NO  VALUES  TO  OBTAIN  VALUES  FOR  MANDATORY  LEVELS 


I JCJ=JC*J28-1 

I IF  ( ABSCOATC JCJFl) ) .GT.SPOMXI  GO  TO  127 

1 ANG=ACR*DATCJCJ ) 

IF  ( I*^  .GT  . ISC  JS  ) . AND.  ABS  (URC  JS)  ) .LT  .SPOMX  ) GO  TO  127 

UR( JS)--OAT<  JCJ  + 1)4SIN( ANG) *CKM 
VRC JS)=-DAT( JCJ41 )*COSf ANG)*CKM 

127  J2=2 

DO  29  K=2.KT 
JK-t  K-1  ) • JT  ♦ JS 
HRJK*HR<  JK )4STDIK) 

DO  28  JH=J2.J17 
JG=JCJ*<  JH-1 ) *2 
JI=JG+1 

IF  CABSCDATl  J I ) > .LE .SPOMX ) GO  TO  128 
IF  ( ABSC OATC J 1-2 )) .GT. SPOMX ) GO  TO  28 
DATC  JG)  = OATI JG-2 ) 

OATC  JI ) = DATC  J 1-2) 

GO  TO  129 

128  IPCJK)=1 

j IF  C ABSC OATC J 1-2 ))  .LE, SPOMX )GO  TO  129 

OATC JG-2 )=D AT C JG ) 

! OAT< JI-2 )=DAT C JI ) 

j 129  IF  CHOZC JH) ♦ESJS.LT ,HRJK ) GO  TO  28 

( IF  I IT.GT . ISC JS ) . AND. IRC JK ) .EQ. 1 ) GO  TO  29 

J 

i 
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RH=  IHRJK-HOZ( JH-1 )-ESJS )/< HOZ( JH )-H02 ( JH- I » I 

UR(  JK  > = (-RH*OAT(  JI  ) «S  IN  ( ACR*0  AT  ( J O ) I 4-OA  T{  J L-2  ) *5 1 N I ACR  *0  AT  IJG- 2 ) ) * 
2<RH-1.0) )*CKM 

VRi JK  > = (-RH«CAT( JI )*COSI ACRAOATI JG>  > ♦OATl J 1-2  I *COS ( ACR*OAT ( JG-2  ) I* 
2(RH-l .0) )*CKM 
GO  TO  29 

28  CONTINUE 
GO  TO  J30 

29  J2=JH 
130  CONTINUE 

IF  ( IT.LT.ISI JSI ) IS(JS»=IT 

30  CONTINUE 
GO  TO  20 

32  CONTINUE 
JSM=0 

DO  ISO  L=1.J1 
IF  ( IS(L>.i.T.lt  )JSM-USM«-1 
150  CONTINUE 

PRINT  l.JSM 

1 FORMAT! 1 OX .37HU  S ARMY  RAWINSONOE  DATA  READ  IN  FOR  t I 2> 9H- STAT 1 ONS 

2 .// J 

RETURN 
END 


I 

5 

'i 
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SUBROUTINE  CHECK  ( J J t YS • XS* VSt US t YL . XL t VN .UN ) 

C 

C THE  SUBROUTINE  CHECKS  THE  OBSERVED  KIND  VALUES  BY  COMPARING  THEM 
C TO  ANALYZED  WIND  VALUES  ITHE  U AND  V OF  THE  INCONSISTENT  DATA  ARE 
C SET  AT  999.9  I. 

C AUG  1977  version:  - MODIFIED  SO  THAT  ANALYZED  WINDS  ARE  BASED  ON  A 
C WEIGHTED  AVERAGING  OF  DATA.  WHICH  IS  MORE  SUITABLE  FUR  USE  WITH  THE 
C PORR  PROGRAM. 

C 

C JJ  = NUMBER  OF  WIND  DATA 

C YS.XS  = LATITUDE  AND  LONGITUDE  OF  WIND  DATA  lOEG) 

C US.VS  = U AND  V COMPONENTS  OF  WIND  DATA 
C YL.XL  = LATITUDE  AND  LONGITUDE  OF  ROWS  AND  COLUMNS 
C UN.VN  = GRID  POINT  U AND  V WIND  COMPONENTS 
C CRT=  CRITICAL  VALUE  USED  IN  TESTING  WIND  DATA 

C W1M=  WEIGHT  GIVEN  TO  A WIND  IN  ANALYZING  A VALUE  AT  ITS  LOCATION 
C AND  WEIGHT  GIVEN  TO  NEAREST  GRID  POINT  VALUE  INOW  CALLED  SIM  AND 

GIM) 

LOGICAL  DEBUG 

DIMENSION  JSTI  26>.0ST(  26).WTS(  25>.OVR(20) 

DIMENSION  XSf  ll.YSI  Il.VSI  1 ) .US(  1) 

DIMENSION  XLC  ll.YLI  1>.UN(  D.VNI  1} 

COMMON/C  IS/  I SI  1000  I 
COMMON/CCHK/  CRT . S I M « G1 M . I WND 

COMMON/CPTS/  KS.Wl .C2 . RMAX . KSSS . I OS * KSW . ALPH 

COMMON/CGO/  M9.N9. I 9 t MB t NB . YB t XB . DO 
EQUIVALENCE  ( USK , USJ ) . I V SK . VS J ) 

EOUI VALENCE  I OVRI 1 > »DNH> . (OVRI2> .OHH) . (OVR(3) .DUH)  .lOVRIA) .OVH) , 

2 (DVR( S) .DTH  ) . (OVR(6) .DXH) . (0VR(7 ) .OYHI . (OVRI 81 .OXYH) . 

3 I OVRI 9 I .OXXhI.IOVRIlO) .OYYH).  (0VR<  1 1 ) .OXHHI. (OVRI 12) .OXUH)  , 

4 IDVRI I 3)  .OXVH) . IDVRI  14) .DYHH) . (OVRI 15) .DYUH) . (DVRI 16) .OYVH) . 

5 (OVRI 17) .OXTH  ) . (DVRI  18 ) .DYTH) 

DATA  DEBUG/. FALSE./ 

OEBUG=.TRUE. 

4 FORMAT(21X. « 1 3H  DELETED  DATA/) 

5 FORMAT  I49H  LATIDEG)  LONIOEG)  UIKTS)  VIKTS)  TEST./) 

6 FORMAT! IX.IOFIO.2) 

PRINT  4 
PRINT  5 

C 

C BASIC  COMPUTATIONS 
C 

JT=JJ 

ACR=3. 1416/180.0 
YDI-1  .O/OD 
XDI=1 .0/DD 
KQ5=KSS5+2 

IF  (IWNO.GT.O)  KQS^KQS^l 
SPM1=ALPH*ALPH 
CR5  = 1 .254CRT 
C 

C ANALYSIS  OF  WINDS  AT  MEASUREMENT  LOCATIONS  (SIMILAR  TO  METHOD  USED 
C IN  SUBROUTINE  MESH). 

C 
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KG=0 
20  JX=0 
I.X=0 
CRH=CR5 

IF  (KG.EQ.2)  CR«=CRT 
80  LX=4.XFI 

IF  (LX-JT)  82.82.100 
82  L=I.X 

IF  IKG.NE.O)  1.=  JSTILX) 

IF  (ABS(US(L> I.GT.IOO.O)  GO  TO  80 
YLM=YS(L) 

XLN=XS<L) 

M=(YU( I >-YLM>*Y0I*1.5 
N=(XLN-XL(1 I )*X0I*1.5 
IF  (M.LT.l) 

IF  (M.GT.M9)  M=M9 
IF  (N.LT.l)  N=1 
IF  (N.GT.N9t  N=N9 
LI=N9*IM-I )FN 
I = ILI-1  >*ICS 
CM=COS( ACRAYLM) 

NOO=0 

K=0 

DO  182  IK=1.18 
182  OVRI1K>=0.0 

IF  (IWND.LE.O)  GO  TO  84 
XSJ^XLI N) 

YSJ=YL<M) 

USJ=UN(LI ) 

VSJ-VNILI ) 

OYS=YSJ-YLM 
OXS=l XSJ-XLN>  *CM 
M^GIH 
GO  TO  89 

84  KsKFl 

IF  (K-KSI  85.85.90 

85  1=141 

IF  INOO-KQS)  86.84.84 

86  J=1SI I ) 

IF  (Jl  84.84.87 

87  IF  I ABS( USI J)  ) .GT.IOO.O)  GO  TO  84 
USJ=US( J ) 

VSJ=VS<J) 

XSJ=XSI  J ) 

YSJ=YS( J» 

DYS=VSJ-YLM 
OXS=<  XSJ-XLN ) *CM 
W=SI  M 

IF  fJ.EQ.L)  GO  TO  389 
0X52=0.0 

DYS2=0VS*0YS4DXS*0XS 

IF  (lOS.LT.l)  GO  TO  88 

OXS1=USK*USK4VSKAVSK40.01 

OXS2=(USK*DYS-VSK4DXS) 

0X52=0X52*0X52/0X51 
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88  W=C2/(0YS2+DXS2*SPHl*C2) 

389  IF  (KG.NE.O)  «3W*NTS(J) 

89  N00=Na04'l 
DNM=ONH+»» 

OUH=OUH4-US  J*W 
DVH=OVH+VSJ*»l 
60  TO  84 

90  IF  (N00.LT.3)  60  TO  80 
IF  (ONH*LE.O.OI  60  TO  80 
DNH=I .O/ONH 
ULL=DUH*DNH 

VLI.=  OVH*DNH 
C 

C CHECK  FOR  FOR  INCONSISTENCY  BETWEEN  ANALYZED  AND  MEASURED  WINDS 
C 

ALTHI=  ULL*ULL-*-VLL*VLL 
ALTH2-  US(L)*US(L)-«-VSfL>*VS(L> 

BLTH=ALTH1 

IF<ALTH2  .6T.  ALTHl)  BLTH=  ALTH2 
PERP1  = (ULL*US<L  ) ♦VLLWVSIL) )/BLTH 
IF  I DEBUG) 

2PRINT  6t  YS(L  ) t XS(L) .USIL) • VS(L) (PERPI (ULL.VLL tCRW .CR5.CRT 
0LS=PERP1 
WTS(L)=1.0 

IF  (OLS.GT.CRS)  60  TO  80 
WTSIL>=O.S 

IF  (OLS.GT.CRW)  GO  TO  80 
IF  (KG.EQ.2)  GO  TO  99 
WTS(L)=SIM 

IF  <JX.GT.25)  GO  TO  100 

JX=JX41 

JST(JX)  =L 

DST(JX)  = PEHPl 

GO  TO  80 

99  PRINT  6.  YS(L  ) t XS(L) *USIL )« VS(L) tPERPl 
US(L)  =-999.9 
VSIL)  =-999.9 
GO  TO  80 
100  CONTINUE 
C 

C ORDER  SUSPECT  DATA  ACCORDING  TO  HIGHEST  PERPl  VALUES  AND  REPEAT  ABOVE 
C 

IF  (JX.LE.O)  GO  TO  300 
JT=JX 

IF  (KC.EQ.OI  GO  TO  202 
IF  (KG.NE.l)  GO  TO  300 
DO  200  J1=1.JT 
J3=J1 

DO  ISO  J2=J1.JT 

IF  (OST(J2)  .LT.  0ST(J3))  J3=J2 

150  CONTINUE 

IST= JST( J1 ) 

XST=OST< J1  ) 

JST( J1 ) = JST(  J3) 

OST( JI )=OST( J3> 
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JST< J3>*IST 
OSTC J3>aXST 
200  CONTINUE 
202  KG^KG*-! 

GO  TO  20 
300  CONTINUE 
RETURN 
END 


r 


SUBROUTINE  AOVECIDT .US . VS t Q 1 • 02 • CM  1 


C 

C THIS  SUBROUTINE  AOVECTS  THE  FIELD  01  WITH  A SPECIFIED  WIND. 

C AUG  1977  version:  - MODIFIED  TO  GIVE  A STABLE  ADVECTION  COMPUTATION 
C UNDER  ALL  INPUT  CONDITIONS. 

C 

C DT  = TIME  STEP  OVER  WHICH  ADVECTION  IS  MADE 

C US. VS  = COMPONENTS  OF  WIND  USED  TO  ADVECT  01  FIELD 

C 01  = FIELD  THAT  IS  ADVECTED 

C 02  = DUMMY  FIELD 

C CM  = COSINES  OF  ROW  LATITUDES 

C 

DIMENSION  OIC  I I .021  I ) .CM< 1 ) 

COMMON/CGO/  M9.N9. I 9 • MB • N8. VB • XB . DO 

XIO=I.O 

YI0=XID*N9 

TT=OT 

UT=US 

VT=VS 

IT=1 

DO  20  IR=1 . I 0 

IF  ( ABS( UT) .LE.l .0 .ANO.ABSI VT ) .LE. I.O)  GO  TO  22 
UT=0.5*UT 
VT=0.5*VT 
TT=0.5*TT 
20  IT=IT^IT 
22  CON  =A8SITT) 

DO  200  IR*1.IT 
1 = 1 

DO  100  M=1.M9 
CMI=1 .O/CMIM) 

DO  100  N=1.N9 
011=01 ( I ) 

UA=ABS(UT>WCMI 
VA=ABS( VT) 

IF  (DT.LT.O.O)  GO  to  80 
IU=N-SIGN(XIO.UT) 

IV=I4-SIGN(YI0.VT) 

GO  TO  90 

80  IU=N4-S1GN<XI0.UT) 

IV=I-SIGN( VID.VT) 

90  01V=0.0 
0IU=0.0 

IF  ( IV.GT.O.ANO.IV.LE.19)  Q1V=01IIV) 

JO=IU*(M-l )*N9 

IF  ( IU.GT.O.ANO.IU.LE.N9)  QIU=01<JU) 

02( 1 )=01 I-l I 01 1-QIU  )*UA*101I-0IV  )*VA)*CUN 
100  1=14-1 

DO  150  I=I.I9 
ISO  01(I)=02(l) 

200  CONTINUE 
RETURN 
END 
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Appendix  B 

LISTING  OF  MODIFIED  AAR  PROGRAM 

PROGRAM  AAR<  I NPOT . OUTPUT » TAPE  1 .T APE 3 . T APE5= I NPU T . 

1 rAPE6=OOTPUTt TAPE? ) 

C 

C PROGRAM  AAR  READS  VARIOUS  INPUT  DATA  DESCRIBED  AS  COMMENTS  IN 

C IN  PROGRAM  and  COMPUTES  A PINAL  FORECAST  CURVE  ( HE  I GHT . V I RTUAL 

C TEMP  IN  DEG.C.  tU  AND  V CUMP.t  FOR  A SPECIFIED  FIRING  TIME  AND 

C FIRING  LOCATION 

C AUG  1977  version:  - CONTAINS  IMPROVEMENTS  SUCH  AS  THE  ADDITION  OF  SUd- 
C ROUTINE  BEST  THAT  PROVIDES  AN  ANALYZEC  GUIDE  PROFILE.  AND  OPTIMIZED 
C PARAMETER  VALUES. 

C 

DOUBLE  PRECISION  XX  , YV , XL . YL . RAO 
INTEGER  ABCOEFI  lAl.GHIJKLI 1 0 ) . MNOPQR ( I 0 ) 

INTEGER  HI<7.7>.TU7.7).U1I7.7».VU7.7) 

DIMENSION  HH(49.4.2).TV(49.4.2).U(49.A.2).VI49.4.2).XM1LS1(21). 
IHHFNAL(4.2) .TVFNAL(4 .2 ) .UF I N AL ( 4 . 2 ) . VF  IN ALC 4 . 2 » . XNOTS 1 < 2 1 ) 
DIMENSION  JSAVIS>.TIM(S).C0IS(21 ) .XNOU(2I ) .XN0V(2I I.XNOTCf 21 ) . 

I XNOT(21 ) .XNOH(2I > . XNOTS2I 21 ) . XMILS2( 21 ) .DIFSEC(21 I 
DIMENSION  UFf  21 >.VFI2I 1 .TF(21 ) .HF(21 ) .PF(21 ) .PFDISI .X(7)  .YI7) 
DIMENSION  UR(12.t2>.VR(12.I2).HR(12.I2).rR(12.12) 

DIMENSION  UZI  144) .VZ< 144 ) .TZ I 144 ). HZ 1 144 ) .PZI 144) 

DIMENSION  IUM(147).XS(2S).YS(25).CS(2S).1M0U2).PR(12.I2) 

DIMENSION  HHTFGIS  ) .STV(S)  . HHTF ( 5 ) . T VTF ( 5 ) . OTF ( 5 ) . VTFI 5 ) 

DIMENSION  US(  12  >.VS( 12) .HS(  12) . TS(  12) .PS( 12) 

DIMENSION  STTFI 12 ) .STPFI 12) .STHF< 12 ) 

COMMON  /FMT/  XLA8 

COMMON/CPTS/  KS.m .C2 . RMAX . KSS5. lOS.KSM.ALPH 

EQUIVALENCE!  I 7 . NX ) . ( J7 . NY ) 

EQUIVALENCE  I UR . UZ  >.( VR. V Z ).( TR . TZ ).{ HR. HZ ).( PR . PZ ) 

DATA  STV/1 1 1 .0.1457.0.3012.0.5574.0.0.0/ 

DATA  IBCK.  ID  IF  . IGMr/0.20.-420/ 

DATA  LI .L4.L12/1  .4.12/ 

DATA  IND. YM AX.  YM IN/21 .300.0.10  00.0/ 

DATA  JSAV/10. 1 .9.10.10/ 

DATA  TIM  /170C* .1830. .2000. .2130. .2300./ 

DATA  PFD/a83.0. 850 .0. 700 .0.500.0.300.0/ 

DATA  IHO/0.31  .59. 90. 120. 151. 181 .212.243.273.304.334/ 

DATA  Ml .C2.KS.KSSS. IOS.ALPH/0.04.0.00S. 10. 10. 1.2.0/ 

DATA  NIL .XN1L/-999.-999.9/ 

DATA  STHF/0.0. 100.0.350.0.750.0. 12 50.0 . I 750 . 0 .2250 .0 . 2750 . 0 . 32 SO . 0 

2 .3750.0.4250.0.4750.0/ 

DATA  STPF/e 74.  0.86 3. 0.837. 0.797. 0.749. 0.70 3. 0.659. 0.61 8. 0.579.0. 

2 542.0.506.0.473.0/ 

DATA  STTF/280. 2 .279.5.277.9.275.3.272.0.268.8. 265.5.262.3.259.0 
2 .255.8.252.5.249.3/ 

2 FORMAT! IX. 19H  NO  DATA  FOR  TFIRE=.F7.1) 

3 FORMAT! IHl I 

4 FORMATI8F10.2  ) 

5 FORMAT! 10X.3F  10  .2 ) 

6 FORMATI2X.12F10.2) 

7 FORMAT! IH  ) 

8 FORMAT!2X.7F10.2) 

10  FORMAT  ! 13A6. A2) 

12  FORMAT!/65H  LATITUDE  LONGITUDE 
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1 ./•2F15.2> 

13  FORMAT  </15X.9HNaRrHlN<ps.020.8»  18H  E ASTI  NSs  tO20«  8) 

30  FORMAT! lX»I2*F7.2.F10*2tF5.2.I5*12XtF8.2.F|0.2iF5»2.I5l 

42  FORMAT(SX.21S.2FI0.2t2I10> 

43  FORMATI5X.3IS1 

44  FORMAT(5X*15.  I10*F10*2) 

45  FORMAT(5X«3I10I 

PRINT  3 
PRINT  1100 
PRINT  3 

THE  NEXT  DATA  SET  IS  INPUT  READ  FROM  DATA  CARDS  IVARIABLES  ARE 
DESCRIBED  IN  COMMENT  CARDS  THAT  FOLLJM 

DEFT=7,5 

READ! St  10)  ABCOEF 
WRITEIBtlO)  ABCOEF 

RE ADI  5 1 30)  IQ.XFIREtXULtXDtNX.TFiRE.YULtYO.NY 
WRITE (6 130)  lO.XFIREtXULtXO.NXtYFIREtYULt YDtNY 
PRINT  1002 
REAOIStlO)  GHIJKL 

READ<5t42)NT.NC«TFCST,TFIREtlOATE. JSTAT 
PRINT  1002 
PRINT  1002 
WRITE(6t 1005) 

PRINT  1002 
READ! St  10)  MNCPQR 
WRITE(6tlO)  MNOPOR 

READ! 5.5) (XS( J)  tYSIJ) .ESIJ)»j3lt JSTAT) 

PRINT  5. IXSI J ).YS( J).ESI J).J-I « JSTAT) 

NEFj=0 

C THE  FCST  DATE  IJDATE)  AND  TIME  (JTIME)  OF  GWC  DATA  IS  NOW  READ  IN 

PRINT  3 

REAO(l)  JOATE. JTIME.NT 
PRINT4StJOATE.  JTIME. NT 
XT  IME=Fl.O ATI  JTIME  ) *0.01 
ITS1=I 

TLATEsTF IRE-  1 .0 
WRITE16. 10 )GHIJKL 

WRITE<6.42)NT.NL,TFCST.TFIRE .1  DATE, JSTAT 

C ABCOEF^HEAOER  LABEL  CARD  FOR  FOLLOWING  DATA 

C lOslNDlCATOR  SHOWING  HOW  FINAL  DATA  IS  TO  BE  PRINTED 

C 1 PRINTS  HARO  CCPY  OF  CURVES  AND  DATA  FOR  CURVES 

C -1  PRINTS  DATA  FOR  CURVES  ONLY 

C XFIRE=LONGITUDE  OF  CENTER  OF  FIRING  PATH  IDEGREES  AND  TENTHS) 

C XUL-LONGI TUBE  OF  UPPER  LEFT  HAND  CORNER  OF  GRID 

C XOsOISTANCE  (DECREES  AND  TENTHS)  BETWEEN  GRID  POINTS  IN  X (LONG.) 

C NXsNUMBER  OF  MINI  GRID  ARRAY  POINTS  IN  X DIRECTION  (LONGITUDE) 

C YF IRE=LATITUDE  OF  CENTER  OF  FIRING  PATH  (DEGREES  AND  TENTHS) 

C VUL=LATITUOE  OF  UPPER  LEFT  HAND  CORNER  OF  GRID 
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C YD  = DISTANCfc  IN  OEG.  AND  TENTHS  OETl»EeN  GRID  POINTS  IN  V DIR.  (LAT» 

C NY=NUMBER  OF  MINI  GRID  ARRAY  POINTS  IN  Y DIRECTION  ILATITUDE) 

C 6MI JKL=HEAOER  LABEL  CARD  FOR  FOLLOWING  DATA 

C NT=NUMBER  OF  I HR  FCST  ARRAY  PERIODS  INCLUDING  TIME  OF  FCST 

C NL=NUMBER  of  LEVELS 

C TFCST=GMT  time  of  FIRST  ARRAY  OF  UPDATED  GwC  FORECAST 

C TFIRE=GMT  TIME  OF  FIRING  (HRS  AND  MINI 

C IDATEsOATE  OF  GWC  FORECAST 

C JSTAT=NUMBER  CF  OBSERVATION  STATIONS 

C XS( J )=LONGITUDE  OF  REPORTING  STATION  IN  DEG  AND  HNDTHS 

C YS( J)=LATITUOE  CF  REPORTING  STATION  IN  OEG  AND  HNDTHS 

C ES( J)=ELEVATION  OF  REPORTING  STATION  IN  METERS 

I Y-INT(  I DATE* 0.0001  ) 

IM=INT(  I DATE A 0.01 )-I YAIOO 
IO=IDATE-IV*l 0000- IM*IOO 
IH=INT(TLATE» 

IMN=( TLATE-IH  )*60 

IMDATE^IMN+bOAI  IHF24*(  IDFlMDf  I M ) *365* I Y-366 I ) - 1 BCK 

REX=0  .0 

Rexi=o.o 

XENO=TFCSTANT-l 

CALL  TIME  ( TF  IRE. TFCST  .TTS. TL  . ITS.  ITL.XENO.HEXl » 

IF  (REX1.EQ.S9)  GO  TO  600 

CALL  CENTERiXFIRE. YFIRE.XUL.YUL .XD.YD.NX.NY.PN.PE.PPS.PW. 

1 I PN, I PE. IPS. I PW. REX) 

IF  (REX. £0.71  ) GO  TO  600 
ITSS=ITS-l 

IFIJDATE.NE.  lOA’-E  (PRINT  1001 
IF( XT IME.NE. TFCST  ) PRINT  1001 
IF( JOATE.NE. IOATE)GO  TO  600 

IF  ( TF IRE  .GT  .XEND. OR. TF  IRE  .LT. TFCST)  PRINT2.TFIRE 

THE  GWC  FCST  CATA  IS  SCANNED  AS  (DUM)  TO  DETERMINE  THE  FCST  DATA 
ON  EACH  SIDE  CF  THE  TIME  OF  FIRING 

IF(  ITSS.LT.  ITSl  ) GO  TO  202 
198  DO  200  JT=ITSl.ITSS 
DO  200  JL-LI.L4 
00  REAO(I)IUM 

IN  THE  211.210  DO  LOOPS.  THE  GWC  FCST  DATA  FOR  LEVELS  LI  TO  L4  IS 
NOW  READ  FOR  THE  FCST  TIME  ON  EACH  SIDE  OF  THE  TIME  OF  FIRING 

202  ITS1=ITSS*3 

PRINT  43. ITSS.ITSl .JT 
DO  211  JT-l,Z 
DO  210  JL^Ll .L4 

READ!  I ) ( (UK  J.  I ) . VI  ( J.  I ) .HI  ( J.  I ) . 1 = 1 . 17)  . J=1  . J7) 

DO  205  1=1.17 
DO  205  J=1.J7 
1H1=INT(H1(J.I)A0.0001> 

Tl<  J.  I ) = I ABS(hl ( J. I )-IHl *1000  0) 

H1(J.I)=IH1 
205  CONTINUE 
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DO  209  JX=1  . I 7 

DO  208  JV=l.J7 

JXV= JX*( JY-1 ) ♦NX 

HH( JXr. JLf JT)=H1( JV, JX)*0.1 

TV( JXYt JLt JT)=T1( JY«JXt«0.1 

U( JXY«JL.JTI-U1( JY. JX)«0.1 

208  Vt  JXY.JL.  JTMVKJY,  JX»*0.l 

209  CONTINUE 

210  CONTINUE 

211  CONTINUE 
JX=IPM 
JY=IPS 

JXY=  JX*-(  JY-  1 ) ♦NX 
JWS=JXY 

jx=ipe 

JXY= JX^I JY-1 )^NX 

JES=JXY 

JY=IPN 

JXY=JX'K  JY-1  I^NX 

JEN=JXY 

JX= IPW 

JXY= JX+l JY-1 I ♦NX 
JWN=JXY 

IN  THE  250  DC  LCOPS.  HEIGHTS*  TEMPERATURES.  U AND  V COMPONENTS  FOR 
THE  VARIOUS  LEVELS  (JL-LI.L4)  ARE  COMPUTED  FOR  THE  LOCATION 
OF  XFIRE  AND  YF  IRE 

DO  250  JT=1.2 
DO  250  JL-L1.L4 

CALL  INTERPIPK.PE.HHI JWS . JL . JT ) . HH ( JES • JL . JT I .XFIRE. 

1 HHXFS ) 

CALL  INTE  RPI  Pl»  , PE  .HH<  JWN  . JL  . JT  ) . HH ( JEN  . JL  . JT  ) . XF  IRE  . 

IHHXFN) 

CALL  INTERPIPPS.PN.HHXFS .HHXFN. YFIRE.HHFIN) 

HHFNALl JL. JT)=HhFIN 

CALL  INTERPtPM.PE.TViJMS.JL.JT) .TV { JES. JL.JT) .XFIRE. 

ITVXFSl 

CALL  INTERPIPW.PE.TVI JWN. JL. JT)  . T V( JEN  . JL . JT I .XFIRE. 

1 TVXFN) 

CALL  INTERPIPPS.PN.TVXFS.TVXFN. YFIRE.TVF  IN) 

TVFNALl JL. JT)=TVFIN 

CALL  INTCRP(Ptl.PE.U(  JWS.  JL.JT)  .U(  JES.  JL.JT)  .XFIRE.  UXFS) 

CALL  INTERPIPW.PE.UI JWN.JL.JTI.Ut JEN. JL.JT). XFIRE. UXFN) 

CALL  INTERPIPPS.PN.UXFS.UXFN.YFIRE .UFIN) 

UF  IN  AH  JL.  JT)=UFIN 

CALL  INTERPi PW.PE. V( JWS. JL. JT ) .V( JES. JL. JT ).XF IRE. VXFS) 

CALL  INTERPIPW. PE. V( JWN. JL. JT ). V( JEN. JL . JT ) ,XF IRE. VXFN) 

CALL  INTERPI PPS.PN.VXFS.VXFN.YFIRE.VFIN) 

VF INAL ( JL . JT  ) = VF  I N 
50  CONTINUE 

IN  THE  260  DO  LOOP.  HEIGHTS.  TEMPERATURES.  U AND  V COMPONENTS  FOR 
THE  VARIOUS  LEVELS  (JL=L1.L4)  AT  LOCATION  XFIRE. VFIRE  ARE  COMPUTED 
FOR  THE  TIME  OF  FIRING  (TFIRE) 
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DO  260  JLsLliL4 

caul  INTERP  ( TTStTU.HMFNAH JL.l  ) f HHFNAL( JL«2 ) t TFIREtHHTF( JL ) J 
CALL  INTERP  ( T T S , TL  . T VFNALI JL 1 1 ) . T VF N AH JL 1 2 ) t TF  IR E . T VTF( JL > > 

CALL  INTERP  ( TT S • TL lUF I N AL( JL . 1 ) • UF 1 N AL ( JL • 2 ) • TF IR E . U TF < JL ) I 
CALL  INTERP  ( TTS t TL . VF IN AL IJL . I » # VF INAU< JL t 2 ) . TF IR E . V TF < JL ) ) 

260  CONTINUE 

00  435  LHT=Ll .L4 

HHTFG<LHT)=HHTFILHT  J4STV(LHT)-ES< I A) 

435  CONTINUE 

liRITE(6.1045)  JOATEtJTIME 
PRINT  1112. UTF 
PRINT  1112. VTF 
PRINT  1112. TVTF 
PRINT  1112.HHTFG 

C SOUNDINGS  OF  THE  AVAILABLE  U.S.  ARMY  RAWINSONOE  STATIONS  (ON  TAPE) 

C ARE  READ  IN  IN  SUBROUTINE  RAMIN 

1=1 

JT=:JSTAT 

CALL  RAM  INI  I . JT  . D I F SEC  .E  S . UZ  . VZ  . HZ  . PZ  . TZ  . I MO.  I MOAT  E.  I DIF.  I GMT  ) 

IN  00  LOOP  311.  THE  BEST  AVAILABLE  STATION  IS  LOCATED  BASED  ON  THE 
TIME  AND  location  OF  THE  SOUNDINGS  AVAILABLE 

DO  311  J= 1 . J T 

XUQO=(XS( J )-XF IRE  )♦!  1 1 .1374C0S<  YF 1 RE*0 . 0 1 745 ) 

YUAO=<VS( J>-YFIRE)*11 1 .137 

CD IS( J )=ABS( OIFSECl J) >*.5*SQRT( XL0D4XL00+YLAD4YLAD ) 

11  CONTINUE 
1A=1 

CM  IN-1  .E»06 
DO  320  J=1 . JT 

IF(CDISIJ) .GT.CMIN. OR. UR(V.J) .LT. -200.0)  GO  TO  320 
CM1N  = C0IS( J ) 

1 A = J 

20  CONTINUE 

PRINT  42. I A 

AT  THIS  POINT.  THE  CURVE  OF  THE  BEST  AVAILABLE  SOUNDING  IS 
PRINTED  AND  THE  CURVES  ARE  MAPPED  FOR  BOTH  THE  BEST  AVAILABLE  AND 
THE  UPDATED  SOUNDING 

KQT  = 2 

3R9  DO  400  L=I.L12 
UF<L)=XNIL 
VF (U )=XNIL 
TF(L)=XNIL 
PF(L)=XNIL 
HF(L)=XNIL 
DO  390  J=1.JT 
US( J)=UR(L. J ) 

VS<  J)=VR<L. J) 

TS( J)=TR(L. J) 
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HS( J)=HP(L, J) 

PS< J)=PH(L,JI 

IF  iOStJ>0LT,-2OOmO)  GO  TO  390 

IF  IPStJ) .LT . 100. 0. OR. PS ( J » .GT. 1000.0)  PStJ)  = STPF(L) 

IF  (TS( J > .LT  .230.0. UR. TS( J) .GT.aoO.O)  T SI J ) = oT TF ( L ) ♦DEF T 
390  CONTINUE 

CALL  8EST(KQT,jT.yS,XS.VS.US,PS.TS.VFIRE.XFIRE. 

2 VFIL > .UEIL » .PE (L » . TFIL)  ) 

HF (L )=STHP  <L) 

393  IF  ( TF (LI  .LT.230. 0 .GR.TF <L  ) .GT.300 .0  ) TF I L I = ST TF I L ) FQEF T 
IF  (PFIL) .LT. 100.0. OR. PFtL) .GT. 1000.0)  PF ( L I = STPF ( L ) 

394  IF  (UF(L)  .GT.-200.0I  GO  TO  396 

IF  (L.EQ.ll  GC  TO  395 
UF(L)=UF(L-l ) 

VF(L)-VF(L-1 ) 

GO  TO  400 

395  UF(L»=200.0 
VF(L)=200.0 
GO  TO  400 

396  UF(L)-UF(L)F200.0 
VF(L)=VF(L)F200.0 

400  CONTINUE 

PFD(  1 )=PF(  I I 
X( I )=XUL 
Y(  1 I =VUL 
DO  410  1=2.7 
XI  I )=XI  I-l  I FXD 
y( I )=y( I- 1 )-yD 
410  CONTINUE 

DO  417  L=L1.L4 
UTFIL )=UTF(L  ) +200.0 
417  VTF(L)=VTF(L)+200.0 

CALL  MOVCURIUF.PF.UTF.PFD.XNDU) 

CALL  FCVCURI VF.PF.VTF.PFO.XNOW) 

DO  425  L=1  .L12 

IF  (UF(L) .LT. -200.0)  GO  TO  425 
UF(L)=UF(L)-200.0 
VF(L)=VF(L)-200.0 
XNOU( L)=XNDU( L ) -200.0 
XNDV(L)=XNDV(L)-200.0 
425  CONTINUE 

DO  427  L=L 1 .L4 
UTFIL )=UTFIL ) -200.0 
427  VTFIL )=VTFIL)-200.0 

CALL  MOVCUR(TF.PF.TVTF.PFO.XNOT) 

CALL  MOVCURI HF .PF.HFTFG.PFO.XNOH) 

IFIIO.LT.l)  GO  TO  450 
XMAX  = 20  . 

XMIN=-20. 

XLAa=6HO-COMP 
PRINT  1010 

CALL  PNTDAT (UF.PF.XNUU.PF . INO , 1 NO . XM AX . XM I N . T M AX , T M I N . bOB ) 
XLAB=6HV-COMP 

PRINT  toil 

CALL  PNTDATI VF .PF.XNOV.PF, I NO. I ND . XM AX . XM I N . T M AX. YM IN . BOB ) 


1 


B-7 


XMAX=3  1 0. 

XMIN=210. 

XUAd=6H  TEMP 
PRINT  1012 

CALL  PNTOATITFtPF .XNCT.PP. INO > I NO . XM AX . XM I N . TH AX . T M I N , BUB) 
XMAX=16000. 

XM  IN  = 0 .0 
XLAa=6HHEIGHT 
PRINT  1013 

CALL  PNTDATI HF .PF  tXKDH  tPF .INOtlNO.XMAX.XMlNtYMAX.YMlN.BOB) 

ABO  CONTINUE 

PRINT  6, (UF(L  } t L = 1 t 12 ) 

PRINT  6. ( VFIL  1 .L=l . 12  ) 

PRINT  6. I TF(L > .L=l  . 12 ) 

PRINT  6.IHF1L > tL=l  . 12) 

PRINT  6, IPFIL  ).L=1 t 12) 

PRINT  7 

PRINT  6.(XN0U(L  l•L=1.12) 

PRINT  6t ( XNOVIL > .L=l . 12) 

PRINT  6. I XNOTIL  ) .L=l  .12 ) 

PRINT  6 . ( XNDHIL ) .L-1 . 12) 

IFIUFI  1)  .EQ.O.O.ANO.VFI  1 ) .EQ.0.0)  GO  TC  AGO 
GO  TO  500 

A60  IF ( TFC 1 ) .EQ.O  .0.  AND.HFI 1 ) .EQ.0.0  ) PRINT  2.TFIRE 
500  CONTINUE 

C WSTM  CONVERSION 

RAD-0.017A532925 

YL=YFIREARAC 

XL=XFIREARAD 

CALL  aSTMI YL.XL.XX.YY) 

PRINT  3 
PRINT  10A6 
PRINT  1002 
PRINT  lOAO 
PRINT  1002 
PRINT  1020 
PRINT  1002 
PRINT  1025 
PRINT  1002 

C IN  DO  LOOP  A7C.  DATA  FOR  BOTH  THE  BEST  AVAILABLE  SOUNDING  AND 

C UPDATED  SOUNDING  ARE  PRINTED  IF.l  FORMAT)  FOR  ZONES.  PRESSURE. 

C TEMPERATURE.  AND  U AND  V COMPONENTS 

DO  A 70  ND=1  . 1 2 

IF  (PFIND ) .GT  .0  .0 ) GO  TO  523 

TF(NO)=XNIL 

UF(NO)=XNIL 

VF(N0)=XN1L 

XNDTI ND)=XNIL 

XNOUINO)  =XNIL 

XNDVINO)  =XNIL 

523  PRINT  1030.  NO . PF < NO ) . TF < ND ) . UF I NO ) . VF ( NO ) . NO . PF ( ND ) . 
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IXNOT(NO) .XNOU(KC) .XNCVtNCI 
470  CONTINUE 

PRINT  1002 

PRINT  12,  YFIBE.XFIRE 
472  CONTINUE 

PRINT  1050 
PRINT  1002 
PRINT  1060 
PRINT  1020 
PRINT  1070 
PRINT  1002 


IN  OC  LOOP  480,  DATA  FOB  BOTH  THE  BEST  AVAILABLE  SOUNDING  AND 
UPDATED  SOUNDING  ARE  PRINTED  IN  COMPUTER  MET  MESSAGE  FORMAT 

DO  480  NO=l ,12 

IFIUFIND) .EQ.C.O.ANO.VFINO) .EO.0.0)  GO  TO  474 
DIR1=ATAN2(-UF(ND) ,- VF I ND))*1 80/3. 14159 
IF  (OIRl.LT.0.0)  DIRI=OIR1+360 
XMILSl (NO  1=0 IRl 4(6400 .0/360.0 1/10.0 
XMPSI=SQRT(UF(NOI*UF(NO)*VF(ND)*VF(NO) ) 

XN0TS1(ND)=1 . 94254 •XMPSl 
GO  TO  475 

474  XMILS1(ND>=:0 .0 
XNOTSl  ( NO  ) = 0 . 0 

75  IF ( XNOU( NO) .EO . 0.0. ANO.XNDV( NO) .EQ.0.0 ) GO  TO  476 
DIR2=AT AN2 I -XNCUI NO ) ,-XNOV(NO ) 1*180/3. 14159 

IF  (DIR2.lt. 0.0)  OIR2-OIR2F360 

XM IL S2 ( NO ) = D I R2 4(64 00. 0/360.0 1/1  0.0 

XMPS2- SORT! XNCUI NO ) ♦XNDU ( NO )♦ XNDV ( ND ) * XNDV ( NO ) ) 

XNOTS2(NO)=l .94254*XMPS2 
GO  TO  478 

76  CONTINUE 
78  N=ND-1 

IXN0T=IFIX( 10.0*XNDT<NO)4.5) 

IPF=IFIX(PF<NC)*.5) 

IXMILl  = IFIX(XMILSl(N0)'*-.5) 

IXN0T1  = IF IX<  XNOTSl ( NO )*.5) 

IXMIL2=IFIX(XMILS2IND  » + .5 ) 

IXNOT2=IF IX( XNUTS2( NO  )4.5) 

ITF=  IF  IX(  1 0 .0  ATFI  NO  ) +0  .5  ) 

IF  (PF(ND) .GT.0.0)  GO  TO  533 

ITF  =NIL 

I XNDT=N  IL 

IXNOTl-NIL 

I XN0T2  = NIL 

IXMIL1=NIL 

IXMIL2=NIL 

33  PRINT  1080.N. IPF, ITF . I XMILI , IXNOTl .N. IPF , IXNDT , IXMIL2 , IXNOT2 
WRITE  (3.1082) 

1 NAR ,N.  IPF  , ITF , I XMlLl  . 1 XNOTl  .N , IPF.  IXNDT . I XMIL2.  IXNOT2 

80  CONTINUE 

PRINT  13.rV,XX 


5 

i 

f 

c 

« 

I 


I 


t' 
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540 

550 

1001 
1002 
1005 
1010 
1011 
1012 
1013 
102  0 
1025 

1030 

1031 
1040 

1045 

1046 
1050 
1 060 

1070 

1060 

1062 

1100 


1110 

1112 

600 


DATA  ) 

OF  graphical 

OF  GRAPHICAL 
or  GRAPHICAL 
OF  GRAPHICAL 


DATA 

DATA 

DATA 

DATA 


fX=U  COMP,  Y=PRESSI/) 
(X=V  COMP,  y=PRESSI/) 
IX=TeMP,  Y=PReSSI/'l 
IX=HeiGMT,  Y=PReSS)/) 


TE 


IF  (ITS2*GT.NT)  GO  TO  550 
DO  540  JT=irS2,NT 
DO  540  JL=Ll,L4 
READ  (1)  lUM 
CONTINUE 
END  FILE  3 

FaRMAT(lX,46H0ATE  TIME  OF  G4C  FCST  DIFFERENT  THAN  DATA  CAROI 
FORMAT! IX, 21 15) 

FQRMAT(17X,17HAHMY  STATION 
FORMAT! IHl ,15X,46HPRINTCUT 
FORMATdHl  ,15X,46HPRINT0UT 
FORMAT! 1 HI  , 1 5 X , 44HPR I NTOUT 
FORMAT! I HI , 15X,46MPRINTOUT 

F0RMAT!4X ,23HeEST  AVAILABLE  5 OUNO I NG ,20X , 1 6HUP DATED  SOUNDING) 
FORMAT! 1 X,70HZ0NE  PRESS  TEMP  UCOMP  VCOMP  ZONE  PRESS 

IMP  UCOMP  VCOMP) 

FORMAT! 2X, I 3, F 7. 1 ,Fe,l,2F7.1,Ill,F7*l,F6.1,2F7.1) 
FORMATI//10X,5HCMIN=,F10.2, 1 OX . 3H I A= , I 5 ) 

FORMAT! 7X«49HPRESS  IN  MBS,  TEMP  IN  OEGIK),  U AND  V COMP  IN  MPS ) 
FORMATIlOX, lOHGttC  DATE  = , I 1 0 , 1 OX , 1 OHGWC  TIME  =,110) 

FORMAT! 16X, 1 SHSTANOARD  FORMAT) 

F0fiMAT!21X,27HCaMPUTER  MET  MESSAGE  FORMAT) 

FORMAT! 1 X ,79hPRESS  IN  MBS,  TEMP  IN  TENTHS  OF  OEGIK).  DIRECTION  IN 
ITENS  OF  MILS,  SPEED  IN  KTS) 

FORMAT! IX .70HZ0NE  PRESS  TEMP  DIR  SPEED  ZONE  PRESS  TE 

IMP  DIR  SPEED  ) 

FORMAT !2X. 13,17,16,217.111,17,16,217) 

F0RMATI2X,11 17) 

F0RMAT1///47X,30HARTILLERY  APPLICATIONS  ROUTINE,  ////// 

2///,  57X,5HLNITS.///,53X,11HSPEED  - MPS,// . 53X , 1 5HD IRECT ION  - OE 
3G.//.  53X, 15HHEIGHT  - ME TE HS . // , 53X , 1 9HTEMPER ATURE  - DEG  K) 

FQRMAT!6X.aHG4C  DATA/) 

FORMAT!4X,5E12.3) 

STOP 

END 
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SUBROUTINE  R «W  I N ( I . J T . T S .E S t UR . VR  .H R . PH . TR . I MO t I D ATE • I O IF . I GMT  ) j 

C THIS  SUBROUTINE  READS  IN  THE  U.S.ARMY  RAWINSUNOE  DATA.  ] 

C AUG  1977  version;  - MODIFIED  TO  SELECT  DATA  TYPES  IT=I  OR  2 ONLY.  AND 
C TO  INTERPOLATE  AND  EXTRAPOLATE  WHEN  WIND  PROFILES  ARE  INCOMPLETE.  ] 

C ! 

DIMENSION  IR(I44)  j 

DIMENSION  UR( I I . VR< I ) .HR( 1 ) .TR(  1 ) ,PR( n .ESI  1 > .TS(  I ),  IMOI 1 I | 

DIMENSION  OAT(SI2).IS(2S).HOZ( 17)  j 

DATA  ACR.CKM, XML/0. 0174533. 0.514B. -999. 9/ 

DATA  J3.J10.J12.J17.J23.J28.J51/3. 10.12*144.23.23.51/ 

DATA  HOZ/0.0 .100.0 .350.0.750.0. 1250.0. 1750.0.2250.0.2750.0.3250.0 
2 .3750*0.4250.0.4750*0.5500.0.6500.0.7500.0.8500.0.9500.0/ 

SPDMX=100.0 

DO  18  L=1 . JT  ! 

TS(L)=XNIL  1 

18  1S(L)=11  I 

DO  19  L=1 . J17  : 

iR(L>=o  ; 

HR(L)=XNIL  ! 

TRtL)=XNIL 

PR<L)=XNIL 

URIL)=XNIL 

19  VRID-XNIL 
IFIL=0 

IF  I I .GT  .1  ) GC  TO  22 
C 

C READS  IN  DATA  SET 
C 

20  READ  (7)  OAT 

IF  (EOF<7  )) 124.21 
124  IFIL^IFILFl 

IF  (IFIL.GT.l)  GO  TO  32 
GO  TO  20 
C 

C LOOP  THROUGH  DATA  SET  IN  RECORD 
C 

21  DO  30  J=1.J10 

22  JC=( J-1 ) • J51 ♦ 1 
C 

C CHECKS  DATA  FOR  D AT E . T I ME . AND  TYPE 
C 

IF  (DAT! JC) .LT.0.0  ) GO  TO  30 
IP1=INT<0AT( JC)*1 .OE-4 ) 

DATJ  = CAT( JC ) - IP  1 • 1 .0E4 
IP2= INTIDAT J*1 .OE-2 ) 

OATJ=DATJ-IP2*1.0E2 
IP3- INTIDATJ ) 

1P4=INT(0AT( JCFl 140*01 ) 

IP5=DAT<JC41 »-IP44100 

JDATE=IP5  + 604(IP4F24  4(  IP3  4 1MDC I P2 1 43654  IP  I -366 ))- I GMT 
IF  I JDATE.LT  . lOATE- lOlF ) GO  TO  30 
IF  ( JOATC.GT . IDATE4 IDIF)  GO  TO  32 
24  JS=INTIDAT( JC42)40. 1 ) 

IF  ( JS.LT.I .OR.JS.GT.JT)  GO  TO  30 
OAT2=OAT( JC42 )-JS410 

1 

I 
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IT=INT(0AT2» 

IF  ( IT.LT.l .QR.1T.GT.2)  GO  TO  30 
esJS=ES( JS» 

TS<  JS)  = JOATE-lt)ATE 

C SELECTS  OUT  PRE SSURE . TEMPERATURE S AND  WINDS  FOR  ARTILLERT  ZONES 
JCJ  = JC«’J28-2 
JX=JC*J3-1 
JKT=l 

DO  28  JH=1 , J23.2 
JK=(JS-l )*J12*JKT 
JJH= JCJ+JH 
HR( JK)=H02( JKT> 

lF(AeS(OAT( JJH+1  ))  .LE.SPOMX)  GO  TO  127 

IF  ( ABS( OAT( J JH-1 ) ) .GT.SPOMX.OR. JH.EQ. n GO  TO  27 

DAT( JJH+1 )=OAT( JJH-1 )S  DAT! J JH ) =0 A T ( J JH-2 ) 

GO  TO  128 

127  1R(JK)=1 

IF  (ABS( DAT! JJH-1 )> .LE.SPDMX .OR. JH.EQ. 1 ) GO  TO  128 
OAT(  JJH-1  1=0  ATI  JJH'*-!  IS  O AT  ( J JH-2 1 =OAT  ( J JH  1 

128  ANG  = ACR*OAT( J JH  1 

IF  ( IT.GT.I S( JS) . AND. IRI JK) .EQ. 1 1 GO  TO  27 
UR(JX)=-DAT{JJH41  )*S1N( ANG) *CKM 
VR( JK)=-DAT( JJH«1 )*COS( ANG) *CKM 

27  IFIDATI JX+JHl  .LT.0.0)  GU  TO  28 

IF  ( I T.GT. IS( JS  I .ANC.TRI JK)  .GT.0.0)  GO  TO  28 
PR( JK)=OAT<  JX4JH) 

TR(  JK)  = DAT(  JX'fJH't-l  M273.16 

28  JKT=JKTF1 

IF  ( I T.LT. IS( JS  1 1 IS<JSI  = IT 
30  CONTINUE 
GO  TO  20 
32  CONTINUE 
RETURN 
END 
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SUbHOUTINE  QESTtKOT.  JJ.YStXS.VS<UStHS*TS«  YL.XL  .VL  .UL.Hi.tTI.) 

C 

C THIS  SUBROUTINE  CCMPUTES  VALUES  FROM  OBSERVED  DATA  AND  AN  INITIAL 
C GUESS  VALUE  (IF  KQT=  2)  BY  A LEAST  SQUARES  FITTING  UF  THE  DATA 

C (SEE  ENOLICH  AND  MANCUSC.  MON  WEA  REV.  1968.  342-350). 

C 

C THIS  SUBROUTINE  WAS  ADDED  TO  THE  AUG  197T  VERSION  OF  THE  AAR  PROGRAM 
C 

C JJ  = NUMBER  UF  WIND  DATA 

C KS  ■=  NUMBER  OF  CLOSEST  DATA  TO  A USED  TO  COMPUTE  ITS  VALUE 

C W1  = WEIGHT  GIVEN  TO  INITIAL  GUESS  VALUE 

C C2  - WEIGHTING  CONSTANT 
C 

C YS.XS  = LATITUDE  AND  LONGITUDE  OF  WIND  DATA  (DEG) 

C US. VS  = U AND  V COMPONENTS  OF  WIND  DATA  (M  SEC- 1 ) 

C UN.VN(KQ=1)  = INITIAL  VALUE  — NOT  USED 

C = FINAL  VALUE  A SMOOTH  ANALYSIS  USED  AS  THE  INITIAL 

C VALUE  IN  The  KQ=2  COMPUTATION 

C UN.VN(KQ=2)  = INITIAL  VALUE  USED  AS  THE  INITIAL  GUESS  VALUE 

C = FINAL  VALUE  USED  AS  THE  FINAL  ANALYSIS 

C HS.HN  = INITIAL  DATA  AND  ANALYZED  VALUE  FOR  AN  ARBITRARY  QUANTITY 
C 

DIMENSION  OVR(20) 

DIMENSION  YS(1).XS(I).VS(I).US(1).HS(I).TS(1) 

COMMON/CPTS/  KS.Wl . C2 . RM AX . KSS5 . lOS.KSW.ALPH 

equivalence  (0VR( 1 > .DNH>  >(DVR(2)«0HH).I0VR(3) .OUH)  . (OVR( 4 ) .DVH) . 

2 (OVR(S ) .DTH ). (0VR(6 ) .0 XH ) « ( D VR ( 7 ) . D YH ) . ( O VR ( 8 ) .OXYH) . 

3 (OVR( 9) .DXXH)  . (DVR( 10)  .DYYH) . (DVRI 1 1 ) .DXHH) . (DVR( 12) .DXUH)  . 

4 (DVR(  1 3) .OXVH  > . (OVR( 14 ) .UYHH) . (OVR(  IS) .DYUH) • ( OVH( 1 6 ) • O V VH ) . 

5 (OVR< 17) .OXTH)  , (DVRI 18)  .DYTH) 

ACR=3. 1416/180. 0 

KS=JJ 

DO  100  KQsl .KCT 
KQ5=KQ-l 4X8 
CM=COS( ACR4YL  ) 

K=0 

IF  (KQ-1)  82.82.83 

82  NOD=0 

DO  182  IK=1.S 

182  OVR(IK)=0.0 
GO  TO  84 

83  DO  183  IK=b.l8 

183  OVR( IK)=0.0 

w=wl 

ONH=W 

NOO^l 

OHH=ML*W 

OUM=UL*W 

OVM=VL*W 

DTH=TL*W 

84  K=K4| 

IF  (K-KS)  85.85.90 


1 

i 


I 
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85  IF  <NCO-KOS>  86.90*90 

86  J~K 

IF  ( ABSIUSIJ ) I.GT .100.0)  GU  TO  84 

XSJ^XStJt 

YSJ*YSt J) 

usj=us( j> 
vsj=vs( j) 

HSJ=HS<  J ) 

TSJ=TS( J) 

OYS=YSJ-yL 

DXS=  (XL-XSJ)4CM 

OYS2=DYS40YS4CXS*OXS 

0XS2=0.5*0YS2 

IF  (lOS.EQ.O)  CC  TO  385 

USK=USJ 

VSK=VSJ 

OXSlsUSKOUSK^VSKAVSK-t-O.O  1 
0XS2=(USK*0YS-VSK*DXS) 
OXS2=0XS24OXS2/0XSl 
385  C2/<DYS24DXS24AUPH+  C2 ) 
N00=N00+l 
ONH  = DNH^Vl( 

TSJ=TSJ*I* 

MSJ=HSJ*W 

USJ-USJ4W 

VSJ=VS J4W 

0HH=DHH4HSJ 

DUM=OUHFUSJ 

OVH=OVH+VSJ 

0TH=0TH>TSJ 

IF  (KQ-I ) 84.84.89 

89  OYH  = OYH+OYS*l» 

0XH=DXHFDXS4»( 

DXYH=DXYH*OXS*OYS*W 
OX  XH=0XXH40XS  40X54 1» 
0YYH=0YYH40YS4DYS4|« 

0XHH=0XHH4HS J40XS 
0YHH=0YHH4HSJ  40 YS 
0XTH=0XTH4TS J4DXS 
0YTH=0YTH4TSJ40YS 
OXOH=OXUH4USJ40XS 
OYUH=OYUM4USJ40YS 
OXVH=OXVH4VSJ40XS 
0YVH=0YVH4VS J4DVS 
GO  TO  84 

90  CONTINUE 

IF  (K0-1>  92.92.94 

92  IF  (ONH)  110.110*93 

93  ONH*l.O/ONH 
HL=0HM40NH 
UL=0UF40NH 
VL=0VH40NM 
TL=DTF40NH 
GO  TO  too 

94  IF  (NOO-31  110.95. 95 
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95  D=DVH*DrH-ONH*OYYH 
E-aDXYH*OYH-DXh*OyVH 
A®OXH#DYM-ONh*DXYH 
B=DXXH*0YH-0XH*0XYH 
BOAE=0*D-A*E 
IF  (BOAE)  97,110,97 
97  a[=1.0/BOAE 

C =DXHH*DYH-DHM*OXYH 
F =OYHH*DVH-OHH*DYYH 
ct=dxth*dvh-oth*oxyh 

FT=OYTH*OYH-OTH*OYYM 
CU=DXUH*DYM-OUH*DXYH 
FU=OYUHAOYH-DUH*OYYH 
cv=oxvh*dyh-cvhadxyf 
FV=OYVH*OYH-OVH*DYYH 
HL-=(a*F-C*E)*EI 
UC=( B*FU-CU*E )*BI 
VL=<  d*FV-CV*E )*BI 
TL=(B*FT-CTFe  )*BI 
100  CONTINUE 
110  CONTINUE 
RETURN 
END 
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Appendix  C 

LISTING  OF  RESULTS 
FOR  INDIVIDUAL  CASES 
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